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By Francois N. Frenkiel 


INTRODUCTION 


In reading the publications on turbulence of different authors, one 
often runs the risk of confusing the various correlation coefficients and 
turbulence spectra. We have made a point of defining, by appropriate con- 
cepts, the differences which exist between these functions . Besides, we 
introduce in the symbols a few new characteristics of turbulence. In the 
first chapter, we study some relations between the correlation coefficients 
and the different turbulence spectra. Certain relations are given by means 
of demonstrations which could be called Intuitive rather than mathematical. 
In this way we demonstrate that the correlation coefficients between the 
simultaneous turbulent velocities at two points are identical, whether 
studied In Lagrange's or in Euler's system. We then consider new spectra 
of turbulence, obtained by study of the simultaneous velocities along a 
straight line of given direction. We determine some relations between 
these spectra and the correlation coefficients. Examining the relation 
between the spectrum, of the turbulence measured at a fixed point and the 
longitudinal-correlation curve given by G. I. Taylor, we find that this 
equation is exact only when the coefficient 

Lx./d^O) 


is very small. 

We find that, in a flow of homogeneous and isotropic turbulence, the 
transverse correlation length is equal to half the longitudinal correla- 
tion length, and we obtain several useful relations with the other char- 
acteristics of turbulence. Next, we introduce nondimens ional parameters 
which greatly simplify the calculations. In the second chapter we view 
a few experimental results and study the method of representing them by 
empirical equations. 


*"Etude statistique de la turbulence. Fonctions spectrales et 
coefficients de correlation." Office National D’liltudes et de Recherches 
Aeronautiques (O.N.E.R.A. ) , Rapport Technique No. 3^> 19*1-8. (The publica- 
tion of this report, completed in 19 * 1 - 2 , was delayed by circumstances 
related to the war.) 
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The following chapters, which form the main part “of .this study, 
provide all the possible representation of correlation and of spectrum 
curves; we compare with them the few measurements already made and intend 
to compare with them the new measurements which we hope to perform in the 
future. We consider therefore this part as a preface to a later report. 

The large number of curves will facilitate the choice of equations which 
will best represent the experimental points, once the test results will be,, 
available. 

For the representation of a correlation carve by a function, one 
might consider the general interpolation methods which permit, for 
instance, giving an approximate representation of- the experimental curve 
by a geometrical or trigonometrical polynomial. Nevertheless, this crude 
method is not favorable because it disregards the typical appearance (bell- 
shaped curve ) of the measured correlation coefficients. Several experimen- 
ters have insisted on the importance of the function R(r) = exp(-Kr) for 
representation of the experimental values. We have learned lately that 
there exist certain theoretical reasons for admitting this particular form. 
In fact. Professor J. Kamp6 de F&riet notified us at hi6 return from his 
recent voyage to the United States that J. L. Doob (ref. 1) has demon- 
strated, with the assumption that the velocity of a particle is a func- 
tion satisfying the following properties, that: 

1. The stochastic process is homogeneous in time , 

2. This process is a Markoff^ - process 

3. The law of probability of u-(t^j,u(t 2 ) is a Gauss law with 
two variables. 

If all these conditions are realized, the correlation coefficient 
between the two velocities of- the same particle at the instants t 
and t + h is of the form R(h) = exp(-K|h|). 

This form of the correlation curve which had already been proposed 
by H. L. Dryden also represents relatively satisfactorily the experimental 
points of the tests made recently by A. A. Kalinske (ref. 2). Thus we 
shall begin our representation with the study of this function. But since 
these conditions are not rigorously realized, it - is of- interest to use the 
functions exp(-Klhl )cp(h) where cp(h) is a geometrical or trigonometrical 
polynomial. 

Nevertheless, since the bell-shaped curve brings to mind the classi- 
cal function of Gauss, we shall study this function, too. The more so. 


"'"That is, if one considers the instants t^"< t 2 < t n _^_ < t n for the 
given values of u(t^_), u(t£) , u(t n _i) , the distribution of u (t n ) 

depends only on u (t n _p). - 
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as one sees - when one represents the second-order moment of the experi- 
mental correlation curve as a function of the first-order moment - that 
(fig. 6) the experimental points are dispersed between the points cor- 
responding to the functions of the form exp(-K|r|) and exp(-Kr^) . We 
studied the functions exp(-Kr n ) with 1 ^ n ^ 2, but after having made 
a- few tests with these curves we stopped losing them because they proved 
to be rather impractical. A fifth chapter is devoted to several applica- 
tions of a series of Hermite polynomials for representing the correlation 
law or the spectrum according to the proposition of J. Kampe de Feriet. 

We represented in n ume rous figures the correlation curves and the 
turbulence spectra which correspond to the selected functions. Several 
applications are made with use of the results of turbulence measurements 
in water and in air. 

Thus one will note the good representation of a transverse correla- 
tion curve measured at the National Bureau of Standards (fig. 19) by the 
equation of the form 


Ry(y) = exp (-Kjy) cos (K^y) 


Two other curves measured in this same laboratory (figs. 25 and 2k), and 
also a longitudinal correlation curve of A. A. Hall (fig. 25) are very 

well represented by the equation R(r) = [* + (1 - Ajcos^rjjexp^KJrl) . 
Other results (figs. 32 and 33) are represented by the func- 
tion R(r) = [l + Arjexp^-KJr \ ) or (figs. 37 to h2) by 

R(r) = Jl + A-jjr + exp | r | ^ . Interesting results are given by 

the function R(r) = A exp^-K-jJrj) + (l - A)exp (-K^ |r | ^. 

The functions derived from Gauss 1 function do not give as interesting 
results j| still, we have been able to represent a few measurements (figs. 7^ 
to 78) quite satisfactorily by the functions of the form 

R(r) = A exp^-K^r 2 ^ + (l - A)exp ), and a correlation curve of 

E. G. Richardson (fig. 60) is rather well represented by a Gauss curve. 

Once the correlation curve is represented by a function, it is easy 
to calculate the different statistical characteristics of turbulence by 
making use of this function. One can thus determine, for a flow of 
homogeneous and isotropic turbulence, the transverse correlation curve 
from the longitudinal correlation curve or vice versa. If the two curves 
have been measured in one and the same flow, it will be easy to verify 
whether an isotropy of turbulence exists in it. We made several comparisons 



4 


MCA TM 1436 


of this type. One -will note particularly the rather remarkable results 
given in figures 33 and 57 where the experimental points are perfectly 
represented by the correlation functions calculated for a flow of 
isotropic turbulence . 

Finally, we represented the spectrum of turbulence measured at the 
National Physical Laboratory (fig. 93) by a curve corresponding to a 
correlation function of the form R^h) = A exp ^ + (l - A)exp(-Kgh2) 

which represents the experiments,! points better than the curve 
R^h) = exp(-K x jhj ) used by H. L. Dryden. 

The possibility of a rapid calculation of the spectral function, once - 
the form of the correlation curve has been determined, is highly important 
since a graphical calculation based on the experimental correlation curve 
is extremely troublesome. - •- — 

Since the work for this report has been done in 1941 and the report 
has been edited in 194-2, we could not take into account” the results found 
by various investigators with which we have only recently become acquainted. 


SYMBOLS 


1. General Symbols for Homogeneous Turbulence 


Oxyz 

reference axes fixed in space. If there exists a 
mean velocity, the axis Ox is parallel to this 
velocity (system of Euler). 

Oxyz 

reference axes fixed with respect to the center of 
a particle. If there exists a mean velocity, the 
axis Ox is parallel to this velocity (system of 
Lagrange ) . 

t 

time 

h 

time interval 

cr 

cyclic frequency (equal to frequency times 2rt) 

t(x,y,z,t) 

value of a scalar quantity at an instant t and at 
a point x,y, z 

*(x,y,z,t) 

mean component of the quantity t taken with respect 
to the time , independently of- the physical proce- 
dure employed 





NACA. TO 


5 


* (x,y,z,t) turbulent component of the quantity ijr 

« ’ 

ty(x,y,z,t) = if(x,y,z,t) + if ' (x,y, z,t ) 


t(x,y,z,t) 



y,z,s)ds 


V (x,y,z,t) = 0 
. r +T 

lim ± / +(x,y,z,s)ds = <p(x,y,z) 

T — »oo ^ J_T 


If cp(x,y,z) = constant, the turbulence is homo- 
geneous. One may assume equality between the time 
space averages. 


lim i f i|r(x,y,z,s)as 

T — >oo J_t 




s,y,z,t)ds 


i r +y 

= lim -±- / i)r(x,s,z,t)ds 
y oo v_/ -y 


r +z 

= lim ^ / i|r(x,y,s,t)ds 

2 . ■ . > 00 w — Z 


[V (x,y,z,t)] 


2 


mean square 



2 




standard deviation 
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t'fxpypZpti) ilr*(x2,y 2 , Z2 ,t 2 ) 

= R^, - 


l/F (xi,yi,zi,ti)] 5 /F ( x 2>y’2^ z 2^ t 2)J 


correlation coefficient 


(in the Euler system) between two fluctuation com- 
ponents of a scalar ^quantity i|r' ( x i , ^i ,z i , ^i) 

and f ( X 2 , '^2 ,Z 2 , ^2)' subscript k depends on 

the relation which exists between these two quan- 
tities. For instance 


sj(x)-4<x). r(x,y,z,t) 


/&• (0,y,z,t)] 2 /[V(x,y,z,t)] 2 



correlation length in Euler's system 


(“) = 


spectrum of turbulence in Euler's system. Represents 
the contribution of the oscillations of cyclic fre- 
quency a) to the mean square 


(x^ZjtjJ 2 

This spectrum may be obtained by harmonic analysis 
of the curve which gives ijf ' (x,y, z,t) as a function 
of one' of the coordinates. The subscript k depends 
on the curve from which the spectrum derives. 




r(X 1 ,Y 1 ,Z 1 ,t 1 ) r(X2,Y 2 ,Z 2 ,t 2 )_ 

/ [r(Xi,Yi,Zi,ti)] 2 j/jy (X 2 ,Y 2 ,Z2,t 2 )J 2 


in the JLagrange system 


between two fluctuation components of a scalar 
quantity (X-^Y^Z-^t-jJ and ^'(X^Y^Zg,^). 

The subscript k depends on the relation existing 
between these two quantities. 
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correlation length in Lagrange’s system 


Tjr 

f KL (<u) 


spectrum of turbulence in Lagrange's system. Repre- 
sents the contribution of the oscillations of 
cyclic frequency cu to the mean square 

gr’(X,Y,Z,t)] 2 


This spectrum may be obtained by harmonic analysis 
of the curve -which gives t|t' (X,Y,Z,t ) as a func- 
tion of one of the coordinates. The subscript k 
depends on the curve from which the spectrum 


derives. 



= 1 


2. Symbols for the Study of the Statistical Relations 
Between the Parallel Velocities 
V 

u,v,v 


instantaneous velocity 
components of the instantaneous velocity 



mean velocity of the flow 


V’ 


turbulent velocity 


u' = u - u 

v’ = V ► 

w’ = w 


components of the turbulent velocity 


A 

r 


direction of a straight line 
measured distance in the direction A 


y|L 



components of the turbulent velocity parallel and 
perpendicular to the direction A 
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correlation coefficient (see table I) 
spectrum of turbulence "'( see table I) 


coordinate at which the correlation coefficients 
between the simultaneous velocities H^(r), R x (x) 

or Ry(y) became zero for the'”first time (fig. 2 ) 


II 

9 

f ° R&(s)ds 

Jo 

II 

9 

r x ° / v 

/ R x (s)ds 

V. 

Jo 

J 

II 

9 

N / 

JT> 

r y o 

/ Ry(s)ds 

Jo 1 


"apparent" correlation length represented by the 
area bounded by the positive part of the cor- 
relation curve measured up to the first point 
for which R = 0 (fig. 2(a)) 


x .% <ap) 


ratio of "apparent" correlation length and true 
correlation length 7 ]' 


[p 2 ] = f s 2 f (s)ds 


s ^y( s )^ s r dispersion of a turbulence spectrum 


't-J 0 
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standard deviation of a spectrum of turbulence 


A length representing the size of the smallest eddies 

responsible for the energy dissipation, in the 
study of the simultaneous turbulent velocities, 
in Euler's system 

At time corresponding to the dimension of the smallest 

eddies responsible for the energy dissipation in 
the study of the turbulent velocity at a fixed 
point as a function of the time in Euler's system 


Nondimensional symbols 


l,l V 

V v* V 

T > T m' T L> T Vm' T L 


dimensionless variables obtained by dividing r, y, 
z, or h by the corresponding correlation length 
(see table I) 



correlation coefficient given as a function of a 
dimensionless variable (see table I) 


a 

cp(fl) 


*o-£ 


quantity representing the cyclic frequency in dimen- 
sionless symbols (see table I) 

turbulence spectrum in dimensionless symbols (see 
table I) 


coordinate for which the correlation coefficients 
B A (r) , ?x( x )j or 5y(y) become zero (fig, 2(b)) 
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p(ap) = r 

t (ap) 


L^ap) 


X 




,(ap) _ y 


(ap) 


PO 



I^ap) J 

(ap) 

r 0 

& 

4* 

II 

(ap) 

x o 


II 

>n 

(ap) 

_ yo 




(ap) 


length, r , x , or y referred to the apparent 
correlation length L^( ap ^ } r (ap) 

(fig. 2 ) 


coordinate for which the correlation coefficients R A . 


R 


■x* 


^ w V.WXJ. J-lwX&JUWO a 

or By given as functions of p( a P), £( a P), 


rj( a P) become zero for the first time 


min 


pH. 


'min 


lp (aP 3 5mln 


J max 


H 

IP 4 ] cp max 


minimum value of the correlation coefficient ? A (p) 

value of p to which the minimum correlation coef- 
ficient corresponds when the correlation curve 
represents B as a function of p 

value of p( a P) to which the minimum correlation 
coefficient corresponds when the correlation curve 

represents R as a function of p(ap) 
maximum value of the spectral function C P A ^ A ) 


value of & A for which the spectral function cp^ ^ 2 ^ 
is maximum 
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dispersion of the spectrum ^(^l)* ® ie sub- 
script "1" corresponds to that of as 

a function of which the spectrum is repre- 
sented, and "2" to the subscript of cpg 



moment of order K of the area bounded by the 
correlation curve R^(p) 



s K Rjs)ds 


moment of order K of the area 
bounded by the "apparent" corre- 
lation curve ^up to the first 

value R^(p) = o) 



moment of order K of the area bounded by the 
spectral curve C P^(^) 





coefficient representing the dimension of the 
smallest eddies responsible for the energy 
dissipation by turbulent viscosity 


<5 

K 

n 

s 


3. Miscellaneous Symbols 

coefficient of viscosity of a fluid 

energy dissipated by turbulent viscosity per 
unit volume 

constant integer 

variable integer 

integration variable 
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c 

m 

A 

B 

a 

P 


e 

log(a) 
exp(a) = e a 


coefficients used for representing the equations of 
the correlation curves 


base of the Napierian .logarithms 

logarithm with the base e 

symbol used for simplifying writing the equations 


erf (a) 


— f exj)(-s2.)ds 
\/jt JO 


numerical function the value of which is to 
be found in the tables 


„ (-l) n (2n) 1 T 

2 n ni L 


2 


4 a 2 + 2 2 

2 .’ kl 


] 


Hermite 

polynomials 


>W»> 
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CHAPTER I 

* 

THEORETICAL CONSIDERATIONS 


1. Correlation Coefficients 

The study of the turbulent motion in a fluid flow requires knowledge 
of the amount of displacement velocities of the macroscopic particles by 
groups of molecules the velocity fluctuations of which present a similar 
appearance. One of the most .important characteristics of a turbulent 
motion is the correlation coefficient which allows numerical definition of 
the s imilar ity of the fluctuations in velocity and better visualization of 
the size of the macroscopic particles. 

The mathematical definition of a correlation coefficient between two 
scalar functions of time and t* 2 is given by the expression 



in which h is a time interval which may be equal to zero when the cor- 
relation between two simultaneous functions ty'pCt), if' 2 (t) is Leing 
studied. The bars represent time averages obtained from measurements made 
at different instants; all other conditions are equal, though. Thus 



lim A- r i|r(x,y,z,s)ds = <p(x,y,z) = C te 
T — »» x J -T 


In order to simplify the study still more, we assume that the time 
averages are equal to the space averages, whether these latter are taken 
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in the entire flow domain or only along an infinite straight line of 
arbitrary direction. One will have for the axes Ox and Oy, in parti- 
cular, the relation 

. n+T 

(x,y,z,t) = lim ~ / (x,y,z,s)ds 

T — >oo J-T _ 


V (s,y,z,t)ds 


(x,s,z,t)ds (2) 


Thus the averages which occur in the definition of the correlation 
coefficient may be taken in still another manner: by repeating the meas- 

urement of ^ 1 ^ 1 2 at different points of the flow, “all other conditions, 

though, being equal. 

The turbulence is considered from the point of view of Euler if one 
studies the velocities at the flow points without concerning oneself about 
the particles which are situated at these points at the instant of the 
measurement . The turbulence is studied from the view point of Lagrange 
if one follows the particles in their motions and notes their velocities 
without concern for the location where they are placed. 

(a) System of Euler .- Assume two particles A and B (fig. 5) the cen- 
ters of gravity of which are, at the instant t, at the points P and Q, 
placed oh a straight line of direction A, at a distance r. The mean 
velocity of the flow is u = U, and the turbulent velocities of these 
particles V'^(t) and V’g(t) are equal to what we call abbreviatedly 

velocities at the points P and Q: V’p(t) and V'^t). After a time 
interval t^ - t the particles will have left the points P and Q, 

after having described trajectories which depend on the joint effect of 
the mean velocity and of the turbulent velocities of the two particles. 

At the instant t^_ two other particles C and D are to be found at 

the points P and Q. The velocities of the particles A and B are 
V^t^) and V'g(ti) 811(1 are different from the velocities at the 

points P and Q: V'p^t^ = V' c (t^) and V'Q^t^ =_ ^'d(^i)* At the 

instant t 2 , the particles C and D, in turn, will also have left the 
points P and Q, and will be replaced, by two other particles: E and F. 




p+X 

= lim 

l 

/ 

X *°o 


/-X 



p+y 

= lim 

JL 

r 

y >oo 

2y, 

'-y 
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In order to obtain in Euler's system a correlation coefficient between 
two s imult aneous turbulent velocities, one studies the velocities at the 
points P and Q without considering the particles situated then at these 
points. This coefficient is given by the expression 


V , p(t)V , Q (t) 

/[ v ’ p(t )] 2 /| j , «(*)] 2 

In order to obtain the averages which occur in this expression, one 
recommences the measurements at different Instants. In the average 

V , p(t)V l Q(t) one will thus take into account the products V t p(t)V , Q(t) , 
V ’p( t l) V 'Q( t l)> V 'p( t 2) V 'Q( t 2)’ 


Since the time averages are supposed to be equal to the space aver- 
ages, one can obtain this correlation coefficient by another method. One 
makes all measurements at the same instant t and calculates the products 
of the turbulent velocities at numerous pairs of points placed on straight 
lines parallel-to A, at distances r, such as V'p(t)V'Q(t) , V'j^CtJV'gCt), 

etc. These products will serve to determine the average V'p(t)V' q(t) . 

One may consider the velocities of the points situated on the same straight 
line as the points P and Q, or on an entirely different straight line. 


L L 

Assume V* „ and V* _ 

A,P A,Q 

bulent velocities in the direction of 


to be the components of simultaneous tur- 

T TM T 


A and V’ 


nents perpendicular to A 


A,P> 

located in the same plane. 


. . _ the compo- 
A,Q * 

One will have two 


correlation coefficients in this Euler system (fig. 1(a) ) . 



The correlation coefficients between the parallel and the perpendi- 
cular turbulent-velocity components at the points P and Q (fig. l(b) ) 
are given by 


u. 


B A ( r > = R A (r) 


tt'A,-p u 'A,q 


/( u 'a,p ) 2 |/( u 'a , q) : 


B>) = 


V *A,P V 'a,Q 


/( 7 'a.p) 2 /( V ' A,qf 
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We call R^(r) abbreviate dly: the correlation_coeff icient in the 

direction A. 

If A coincides, in the relations (3), with the x axis, parallel 
to the direction of the mean velocity (fig. 1(c)), one will write 

R^(r) = R^(x) = R x (x) and R^(r) = R^(x) (3* ) 

R x (x) is the coefficient of longitudinal correlation between the 
longitudinal turbulent velocities, and we shall call it abbreviatedly: 

longitudinal-correlation coefficient. , B^(x) is the longitudinal- 
correlation coefficient between the transverse turbulent velocities. 

When A coincides with the y-axis perpendicular to the direction of 
the mean velocity, one has 

R J r ) = R y(y) = By(y) R ^(r) = R J(y) (3") 

where Ry(y) is the transverse-correlation coefficient between the longi- 
tudinal turbulent velocities - we shall call, it transverse-correlation 
coefficient - and R^(y) represents the coefficient of transverse cor- 
relation between the transverse velocities. 

One may define, in Ruler's system, a correlatiori“coeff icient between 
the turbulent velocities at the same point but at two instants t 
and t + h. This coefficient is given by the expression 


V , p(t)V'p(t + h) 



One obtains the avera ges by varying the initial instant t. In the 
average V , p(t)V'p(t + h) one will thus take into account the products 

V'p(t)V'p(t + h) , V'p(ti)V' p (t 1 + h). Applying equation (2), one may 

form the means by measuring the velocities at two instants t and t .+ h 
only, but taking the measurements for an infinite number of points. In 
this case the products V , p(t)V'p(t + h), V'Q(t)V'Q(t + h), etc. would 

occur in the calculation of V'p(t)V , p(t + h) . 
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Studying the correlation between the components of the turbulent 
velocities parallel to the direction of the mean velocity, we shall have 


Rt(h) = 


u* (x,y,z,t) u» (x,y,z,t + h) 
(x,y,z,t)] % (x,y,z,t + h)] 




which we shall call abbreviatedly: correlation coefficient at a fixed 

point. For the components perpendicular to the direction of the mean 

velocity, one will have the coefficient R^(h). 

The correlation coefficient (4) has been studied by taking a point 
fixed with respect to space. One may also determine a correlation coef- 
ficient for a point being displaced with the mean velocity of the flow, 
and one will have 


RJ 11 ) - 


U* (x,y,z,t)u* (x + Uh,y,z,t + h 


/&* (x,y,z,t)j (x + Uh,y,z,t + h)] 


(5) 


This coefficient, as also the correlation coefficient R^(h) between 
the transverse components, appertains to a pseudo-Eulerian system. 


(b) System of Lagrange .- Let us now study the correlation coefficients 
of Lagrange’s system. Assume QAjAgAj to be the coordinate axes the 

origin of which is constantly located at the center of gravity of the 
particle A and which are situated so that the axis QAj_ is always 

parallel to the direction A. The correlation coefficient of the Lagrange 
system between the simultaneous turbulent velocities of two particles 
placed on a straight line in the direction A, at a distance r, will be 


V’ A (0,0,0,t)V’ (or,pr, 7 r,t) 


where a, p, 7 are the direction cosines of the straight line A. 

In the average V'^COjOjOjtJVKo^prjTTjt) products like V , A (t)V l B (t) , 

V* ^(t^V^^t^ , V'A^iyv^H^i), etc. will occur. With the assumption 

that the averages with respect to time are equal to the averages one would 
obtain If one made an infinite number of measurements at the same instant t. 


I (y’ A ( 0 , 0 , 0 ,t )] 2 W |v’(ar,pr, 7 r,t)] 
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which amounts t o supposing the exis tence of an ergodic principle, one may 

obtain V , ^(0,0,0,t)V l (<xr,(3r,7r,t) in another manner: one takes the 

averages of- the products of turbulent velocities of pairs of particles 
placed, at the instant t, on straight lines of direction A and at dis- 
tances r so that, for instance, V’^(t )V' B (t) , V'j^tJV' ^(t ), etc. One 

now has the equalities V'^(t)V , -g(t ) = V , p(t)V , Q(t) and 

= V'p(t)V , g(t) . Consequently, the correlation coefficient 

between the simultaneous turbulent velocities of the particles placed on 
a straight line in the direction A and at distances r is obtained in 
the same manner as the correlation coefficient between the simultaneous 
turbulent velocities at the points placed oh a straight line of the same 
direction and situated at the same distance. When the turbulence is 
homogeneous, the correlation coefficients between the simultaneous turbu- 
lent velocities for a given direction and distance are the same, whether 
the study is made in Euler’s or in Lagrange’s system. 


R A , E (r) = H A ^ L (r) _ (6) 

The correlation coefficient between the turbulent velocities of the 
same particle is given by the expression 


v *A( t ) v 'A( t + h) 

+ h ^ 2 

where the average is taken by varying the initial instant t. Therefore, 
the products V'A(t)V’^(t + h), Y'^ t]_ V'^t^ + h), among others, will 

occur in V’^(t)V'^(t + h) . This correlation coefficient may be obtained 
also by studying the turbulent velocities of an infinite number of parti- 
cles at two instants t and t + h only. The average ^’^(tjv'^t + h) 
will then be calculated with the products V'^(t)V'^(t + h) , 

+ h), etc. Thus, if one has the equality 

V’ A (t)V’ A (t + h) = V'p(t)V' p (t + h), one will, on the other hand, have 
in general the inequality V'p(t)V'p(t + h) f V'Q(t)V’Q(t +-h). Conse- 
quently, the correlation coefficient between the turbulent velocities 
at a fixed point (system of Euler) and the correlation coefficient between 
the turbulent velocities of the same particle (system of Lagrange) for the 
same time interval are not necessarily equal when the .turbulence is 
homogeneous . 
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The correlation coefficient "between the components (parallel to the 
direction of the mean velocity) of the turbulent velocities of the same 
particle A , at two instants t and t + h, will be equal to 


u u* A (t)u' A (t + h) 

V h) - B tL (h) - 


^[u* A (t)] 2 | ^|u* A (t + hff 


(T) 


and ve shall call it abbreviatedly correlation coefficient in the Lagrange 
system. The analogous correlation coefficient between the components per- 
pendicular to the direction of the mean velocity will be designated by 

Rti( h ) • 

(c) Homogeneous and isotropic turbulence .- In the preceding section 
we have assumed that the turbulence is homogeneous ^ that is, that a trans- 
lation of the axes does not produce a change in the value of the averages. 
If one assumes furthermore that the turbulence is isotropic, a rotation 
of the axes will not have any effect, either, on the value of the averages. 
Hence there result the equalities 


v v 

R x = *y *y 555 ^ 


( 8 ) 


2. Spectra of Turbulence 

The turbulent energy of a fluid medium may be considered as the sum 
of the energy of simple harmonic vibrations of different frequencies. The 
character of the turbulence will be completely defined if one knows the 
total turbulent energy and the proportion of energy corresponding to every 
frequency, that is to say, the spectrum of turbulence . If, in the study 
of a flow, the observer displaces himself with the mean velocity of the 
flow, it is easier to study the spectrum as the energy distribution as a 
function of the wave length. Every turbulent medium may have spectra 
peculiar to it ■which will not be similar to the spectra of another medium 
except for the particular cases where for instance the causes which pro- 
duce the turbulence are similar. 

One can see an analogy between the spectrum of turbulence and the 
spectrum of light. It must however be noted that, in the case of turbu- 
lence, the spectrum may be subject to a transformation which does, not 
directly depend on external causes but which is due to the production of 
the small eddies initiated by the large ones. The spontaneous decrease 
in intensity of the longitudinal turbulence downstream of a grid is one 
of the consequences of this phenomenon. 
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In order to explain what the different spectra which we shall intro- 
duce in this study are representing, we shall assume availability of a 
measuring apparatus capable of recording the components of the turbulent 
velocity at every instant and simultaneously for several particles. Let 
one recording represent the values of the longitudinal components of the 
instantaneous velocities of the particles which are at a given instant 
situated on a straight line parallel to the direction of the mean velocity. 
The longitudinal spectrum of turbulence may’ be determined by ™.ld ng a 
harmonic analysis of such a recording. Instead of studying the distribu- 
tion of the simultaneous velocities along an infinite straight line, one 
can take several recordings for straight lines of a length which is suf- 
ficient to make the correlation corresponding to this length negligible. 
For each recording, one makes a harmonic analysis and then determines the 
mean distribution of the longitudinal turbulent energy (proportional to 

u' 2 ) which corresponds to this spectrum. Thus, the longitudinal spectrum 
of turbulence f x (cu) will be obtained. By recording the longitudinal com- 
ponents along the straight lines located in a plane perpendicular to the 
direction of the mean velocity, one obtains the transverse spectrum of 
turbulence f(oi) . 

If one makes the harmonic analysis of a recording which represents 
the longitudinal component' of the turbulent velocity at a fixed point as 
a function of the time, one obtains the spectrum of G. I. Taylor f ^ (to) . 

These three spectra correspond to Euler's viewpoint. 

The recording which gives the longitudinal component of the turbulent 
velocity at a point which is displaced with the mean velocity of the flow, 
leads to a spectrum of a pseudo-Eulerian system which we call "spectrum 
following the mean motion f m (co) 

Finally, one obtains a spectrum in the Lagrange system by taking as 
a basis a recording which gives the fluctuations of the longitudinal com- 
ponent of the turbulent velocity of a particle. This spectrum will be 
called spectrum of J. Kampl de F4riet f^ L (o 3 ) . 

In an analogous manner, one may obtain spectra which give the dis- 
tribution of the transverse turbulent energy (proportional to v ' ^ ) by 
starting from the recordings which represent the fluctuations of the 
transverse component of the turbulent velocities. Thus, one will have 
the spectra which will be designated by — - 

*>> *>> *>) *>) < T 

x y t m tL 

To our knowledge, the only direct measurement s of spectra made so 
fax are measurement s of spectra of G. I. Taylor. 
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If the turbulence is homogeneous and isotropic, there correspond to 
the equations ( 8 ) the relations 


f x (o)) = fy(co) 


f y (co) = f x (o)) 


( 9 ) 


5 . Relations Between the Spectra of Turbulence and 
the Correlation Coefficients 

In reference G. I. Taylor studies the relation existing between 
the spectrum of turbulence measured at a fixed point and the correlation 
between the simultaneous longitudinal velocities at two points situated 
on a straight line parallel to the direction of the mean velocity (longi- 
tudinal correlation) . 

Assuming that the turbulent velocity is very small relative to the 
mean velocity (ref. 7 - equation (l)) , he finds that the spectrum and the 
correlation are determined from one another by Fourier transform, according 
to the equations 





cos ^-~^R X ( s ) ds 


( 10 ) 


R x (x) 


r 00 

= J cos(SE)f t (s)ds 


(H) 


Taylor* s calculation may be applied for determining the relation 
betveen the longitudinal spectrum and the longitudinal correlation^ but 
in this case it is not necessary to make a hypothesis concerning the 


magnitude of the turbulent velocity, 
to those of Taylor 



We obtain thus equations analogous 
cos^^R x (s)ds (12) 

(vK (s > ds 


(13) 
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Likewise, one may write for the transverse spectrum and the transverse 
correlation the equations 


_ 00 

■ w/ 0 =° e (ir)y 8 > ,is 

(14) 

O °° 

Ry(y) = cos(^f y (s)ds 

( 15 ) 


Study of the turbulent velocities at a fixed point in space will lead 
to the equations 


f-fc(oj) - — cos(o>s)Rt(s)ds 

* J 0 

( 16 ) 

r“ 

R-fc(h) = / cos(sh)f t (s)ds 

Jo 

(17) 


which give the relation between the spectrum of G. I. Taylor and the cor- 
relation at the fixed point. - 

Comparing the experimental spectrum f-fc (oo) with the one obtained by 
application of equation (10) to the correlation curve R x (x) (ref. 7 - 
fig. 1), and inversely, making the comparison between the correlation curve 
and the one given by equation (ll) starting from the spectrum (ref. ^ - 
fig. 2), Taylor has found that the points calculated by Fourier transform 
are very satisfactorily located with respect to the curves. Applying (10) 
and (12), we find 

ft(co) = f x (a>) (18) 


and (ll) and ( 17 ) give 


R t (h) = R x (hU) ( 19 ) 

This comparison does not verify in general the exactness of the equa- 
tions (10) and (ll) as one might be tempted to believe, but only the fact 
that in this particular case the spectrum measured at a fixed point 
the longitudinal spectrum are identical. 
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A. A. Kalinske and E. R. van Driest have made measurements of cor- 
relation coefficients between the transverse turbulent velocities in 
water (ref. 14) . Comparing the correlation curve between the simultaneous 

velocities R v (hu) with the curve for a point fixed in space R v (h) , they 

X *fc 

find that the tvo curves diverge more and more "when h (or x) increases 

(fig. 4). 

This is easily understandable. Thus we consider a turbulent flow of 
the mean velocity U. Assume two points P and Q located on a straight 
line parallel to the mean motion, separated by a distance x, with Q 
downstream of P. The correlation coefficient between the parallel com- 
ponents of simultaneous turbulent velocities at these two points is R x (x) . 
We ass ume now that P is displaced with the mean flow. After a time 
interval h = ^j, the point P will be at Q. Let R m (h) be the correla- 
tion coefficient between the parallel components of the turbulent velocity 
at a point which is displaced with the mean velocity. The correlation at 
a point fixed in space (here the point Q,) depends on these two correla- 
tions. If R m (h) = 1, the two coefficients R-^h) and R x (hU) axe equal. 

Since the correlation R m (h) diminishes in general when h increases, 
the coefficients R^(h) and R x (hU) will differ the more, the larger h 
will be. 

The correlation curve R m (h) obtained by E. G. Richardson (ref. l 6 ) 

following the mean flow of the water does not correspond, either, to the 
curve R x (x) (fig. 5) • Although for this experiment no correlation curve 

at a fixed point is available, one may say with certainty that - like in 
the tests of KalinBke and van Driest - an important difference will exist 
between Rj-(h) a* 113 - R x (hU) because R m (h) decreases appreciably even 

for small values of h. 

For these two series of experiments made in water (refs. 15 and 1 6 ), 
one will thus have the inequalities 

R^(hU) £ R^(h) R x (hU) j- R t (h) 
unless h is very small. 

Consequently, the equations (l 8 ) and (19) which are verified for 
the tests studied by G. I. Taylor (ref. 7 ) will not be verified for these 
two test series. The equations (10) and (ll) are therefore not exact in 
these cases. 
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One may assume that the two curves R x (hU) and Rt(h) are identi- 
cal when the correlation following the mean flow R m (h) remains very 
large tip to. values of h = £ to which corresponds a negligible correla- 
tion between the simultaneous velocities. This amounts to imposing three 
conditions : 


(a) RjnCh) must decrease very slowly, that is, the radius of curva- 
ture at the peak of the curve representing this coefficient must be very 

d2 Rm(q) 


large, and the second derivative at the origin 
absolute value. 


must have a small 


dh £ 


(b) When x increases, R x (x) must become rapidly negligible so 

that to small values of this coefficient still correspond coefficients 
R m (h) of the order of unity. For the correlation curves, the negative 

part of which is negligible, this may be expressed by requiring that the 

p oo — 

correlation length L^. = / R x (s)ds be very small. 

J 0 


(c) Finally, it is necessary that small values of h correspond 
to large values of x or, in other words, that the mean velocity is 
high. 


Thus one will have the general condition that the coefficient 


LjJd^O) 


( 20 ) 


must be- very small. The value of this coefficient could serve as a 
criterion of whether the equations (10) and (ll) may be legitimately 
employed. - 

In the case of the tests studied in the" references 5 and 12, no 
curves R m (h) are at disposal, and the value of (20) cannot be calcu- 
lated. In contrast, it is possible to compare the values of the ratio 

which in the tests studied by Taylor varies between 0.022 and 0.005 second 
and is in Kalins ke' s and van Driest' s experiments equal to 0.118 second. 


Because the number of tests made is insufficient, one can only say, 
in an entirely arbitrary manner, that Taylor's equations may be considered 

exact when the ratio ^ is of the order of 0.005 second, and one must - 

expect that they will not be verified if this ratio is much larger. 
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Studying the turbulence from Lagrange's viewpoint, J. Kampe / de Feriet 
has obtained (ref. 3 ) the relations 


2 f” 

f tL (a>) = -J cos(cos)R tL (s)ds 


( 21 ) 


= f eos(sh)f tL (s)ds 

Jo 


( 22 ) 


which determine the correlation between the longitudinal turbulent Veloc- 
ities of the same particle at two instants t and t + h, with the initial 
instant t being variable. 


4. Equations Capable of Representing the Laws of Correlation 
and the Spectral Functions 

It is not possible to determine a general law for the spectral func- 
tion. However, it is convenient to represent the spectral function in 
the form of a simple equation which evidently can be only an empirical 
expression. Instead of representing the spectrum by an equation, one 
may represent the correlation curve which can be more easily determined 
experimentally than the spectrum. 

For representing the correlation law, the National Bureau of Standards 
has suggested (ref. 11) use of the relation 

Ry(y) = exp^- (23) 

J. Kampe de Feriet studies (ref. 3) & certain number of spectral functions 
which lead, for the correlation law among others, to the law of Gauss 

H^h) JS0 (24) 

and to a lav similar to the one given in the reference 11. 

^(h) = 


( 25 ) 
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These laws are such that the correlation coefficient remains con- 
stantly positive. However, several experimenters obtained negative coef- 
ficients by measuring the correlation between the turbulent velocities 
at two points relatively far apart. This form of the correlation law is 
an indication of important changes in the general character of the 
phenomena due to turbulence, and more particularly in the turbulent dif- 
fusion. Therefore it will be useful to dispose of the correlation curves 
which yield negative correlation coefficients. 

The correlation function Rt(h) and the spectral function f^(co) 
are connected by the equation 


Rt(h) = f cos(sh)f t (s)ds (17) 

J 0 


But there exists a very important difference between these two functions. 
Actually, it is necessary and sufficient that a function, in order to be 
capable of representing a spectrum, satisfies the conditions 

ft(«) > 0 (26) 



(27) 


In contrast, it is difficult to recognize whether R-fc(h) is a cor- 
relation function, that is to say, whether it is the Fourier transform 
of a positive function. These functions form the class of "positive 
definite" functions which have been the object of numerous studies, in 
particular of the fundamental report of S. Bochner (ref. 12). Since the 
criterion formed by the necessary and sufficient conditions which 
S. Bochner presents is not easily manageable, we are, in general, con- 
tent with an examination whether the functions we intend to select 
satisfy the following four conditions which are necessary but not 
sufficient 


-1 Rfc(h) < + 1 

(28) 

r Rt(s)ds > 0 

Jo 

(29) 


lim R t (h) = 1 
h — - > 0 


( 30 ) 
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Urn Rt(h) = 0 (31) 

h — >00 

One may, moreover, require that the first derivative of the correla- 
tion function be zero for h = 0 ■which - with equation (30) taken into 
account - will give 

< 0 


This inequality is not verified for the functions studied in 
chapter III. 

After having assumed a function for representing the law of correla- 
tion, one can calculate the spectral function by making a Fourier transform. 

In order to determine directly the equation which can represent the 
spectral function, one will, beside the equations ( 2 6 ) and ( 27 ), also set 
up the conditions 

__ 00 

lim f t (o>) = | j R t (s)ds (33) 

cd — »0 71 Jo 

lim f t (co) =0 (3*0 

CD 

For each of the functions suggested for representing a correlation 
curve or a spectral curve, the coefficients will be limited in such a 
manner that these conditions will be verified. 

All conditions given for the correlation law R^(h) and for the 
spectrum f-^ (cd) can be applied also to the other spectra and correlation 
laws. In this case, one must replace, in the equations (2 6 ) to (3*0, the 

expressions h, R-t(h), f t (co) by R x (x), Uf x (cD) or by R y (y), 

XJfy(o)) if the correlation between the simultaneous longitudinal turbu- 
lent velocities is studied, by h, R^h) , f m (cb) if the turbulence 

along the mean flow is studied, and by h, R^l, if th e turbu- 

lence is examined from the viewpoint of Lagrange, by studying the fluc- 
tuations of velocity following the particles. 



d%(0) 

dh 2 
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5. Theoretical Equations in a Flow of Homogeneous 


and Isotropic Turbulence 

Ei. Karmah has demonstrated (ref. 10 ) that one has in a flow of 
homogeneous and isotropic turbulence 


B g(r> - #»> ♦ i r 2%) 

Applying the equations (3') and (5"), one finds 


( 35 ) 


%(y) - R x (y) + £ y (36) 

where R x (y) represents the correlation law R x (x) in which the x 

has been replaced by y. This equation gives the relation between the 
longitudinal and the transverse correlation for the longitudinal turbu- 
lent velocities. 


Application of the equations (3 ' ) and (3") yields for the correla- 
tions between the transverse turbulent velocities the relation 


„v, \ v / N 1 dR v (x) 

R (x) = R (x) + - X — — 
x N y w 2 dx 


(36*) 


Integrating the differential equation with the second member ( 36 ) 
and calculating the integration constant -by means of - the condition ( 31 ), 
one finds the longitudinal correlation as a function of the transverse 
correlation _ 


R x(x) = sRy(s)ds (37) 

The tensor equation of Harman (ref. 10, equation (l)) permits cal- 
culation of the correlation R A between the longitudinal turb ulen t 
velocities at two points placed on a straight line of arbitrary direc- 
tion A, as functions of the correlations R x and Ry. One obtains 

r a(|/x 2 + y 2 ) = - g x2 Rx (/x 2 + y 2 ) + R y (/x 2 + y 2 ) ( 38 ) 

x* + y* XT' +-y* 

One may represent- R A as a function of the longitudinal correlation 

or of the transverse correlation alone, by applying the equations ( 36 ) 
or ( 37 ). 
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Equation (13) gives 

3 ®x( x ) = J^q cos (■xr)^x( s ) < ^( sx ) = J^q cos (u)^^)^ 8 

and, taking the condition ( 30 ) into account, one obtains 


lim xEl^Cx) = 0 

x — 

As a result, one finds by integrating the two members of equation ( 36 ) 
the very important relation 2 

W * \ ** ( 39 ) 

In a flow of homogeneous and isotropic turbulence the length of the 
transverse correlation equals half the length of the longitudinal 
correlation. 

In the same manner one obtains by integrating the equation ( 36 ' ) 


v 

L 

x 


= Jl v 
2 y 


( 39 ') 


Replacing Ry(s) in equation (lb) by its value given in ( 36 ), and 
taking (12) into account, one finds the equation 


fy(») = \ - 



df x (o)) 

das 


(bo) 


which gives the transverse spectrum as a function of the longitudinal 
spectrum for the longitudinal turbulent velocities, and applying now the 
equations ( 9 ), one obtains the relation between the spectra which give 
the distribution of the transverse turbulent energy 


f x ( CD) = | fy(cD) - 



cLCyfa) 

dm 


(,b O') 


2 ThlG report was ready for publication when we learned from the 
memorandum of K. Wieghardt in the Luftfahrtf orschung of February 28, 19*1-1, 
that this relation had already been demonstrated. 
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Integrating the differential equation with the second member (40) 
and determining the constant by the equation (33) > one finds the longi- 
tudinal spectrum of turbulence, the expression 

f x (o>) - au P ds (M) 

J CO 

As G. I. Taylor has shown (ref. 6), the mean energy dissipated by 
turbulent viscosity in the unit volume of a flow of homogeneous and 
isotropic turbulence is equal to 


® = 15 


|iU ,g 

X 2 


where A is a length representing the dimension of the smallest eddies 
responsible for the dissipation. This length is linked to R (y) by 
the relation ^ 


-4 = lim 
A y >0 


1 



When a parabola -which passes through its peak is superimposed on 
the transverse-correlation curve, A is the abscissa of the point of 
intersection of the parabola with the y-axis. 

One may also write (ref, 10) the relation 

_1_ = d 2 R x (0) i dgRy(O) 

dx^ ^ dy^ 



which shows that A is equal to the radius of curvature at the peak of 
the longitudinal-correlation curve. In order to obtain a finite dissi- 
pated energy, the second derivative must be zero at the peak of the cor- 
relation curve. 


We call the "dispersion of a curve with respect to an axis" the 
quotient of the second-order moment of the area bounded by that curve 
and the area itself. The standard deviation is equal to the square root 
of the dispersion. Since the area bounded by the spectral curve is equal 
to unity, one will have for the dispersion ;of the longitudinal turbulence 
spectrum the expression 
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and since, on the other hand (ref. 3, p. 172 ) 



a^xCo) 

dx^ 


one finds, applying equation (42) , that the dimension A is equal to 
the ratio of the mean velocity and the standard deviation of the longi- 
tudinal spectrum 


A = 



(43) 


The dimension of the smallest eddies responsible for the dissipation, 
referred to the longitudinal-correlation length, may be given as a func- 
tion of the spectrum alone by the relation 



1 



f x (cu) 


(43*) 


Making the same 
the dimension ~K as 


calculation for the transverse spectrum^ one finds 
a function of the standard deviation of that spectrum 



The equations (43) and (44) give the relation 


(W) 




(45) 


The dispersion of the transverse spectrum is twice the dispersion 
of the longitudinal spectrum when the turbulence is homogeneous and 
isotropic. 


When one studies the correlation between the longitudinal velocities 
at a fixed point, one can define a quantity which is analogous to the 
dimension A and is given by 



c^MO) 


dh c 


( 46 ) 
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The time is equal to the inverse standard deviation of the 

spectrum of G. I. Taylor 



and referring this time to the correlation time L^, one finds 


(* 7 ) 


= 2 1 

11* ft M 

t o >— >0 


( 47 *) 


6. Nondimensional Notations 

In order to simplify the calculations and the notations,, we employ, 
in what follows, dimensionless coefficients. We put in Euler's system 


P = 


IA 


*-£ 




T = 


_ h 


and 


cuL a 

A = ™ 

A u 


CliLy 

CLr m —A 

^ u 




% = 


The correlation coefficients between the simultaneous turbulent 
velocities are 


g A (p) = R A (r) R x (l) = R x (x) 
and the spectral functions 


, » Uf A (oi) , v Uf 




S y (Tl) = Ry(y) R t (r) = R t (h) 


MM • ^ 




%(%) - 
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The dime nsion of the smallest eddies responsible for the dissipation 
can be referred to the longitudinal-correlation length or to the 
transverse-correlation length 


_A_ 




In Lagrange's system one puts 


= 

Lfc 


T L 


h 

LtL 


and 


^tL ■ 


The correlation coefficient between the turbulent velocities of the 
same particle at two Instants t and t + h Is 

5 tL M - B tL (h > 

and the spectrum of J, Kampe de Feriet 


^tL (^tL ) 


f t L (cD) 

LtL 


The conditions which must be imposed so that the function represents 
a correlation law, are written with these notations 


-i c r a (p) + i 
P 00 

m 

/ R A (s)ds = 1 

J 0 

(^9) 

lim R A (p) = 1 

P — >0 

(50) 

lim R A (p) = 0 
P — >o 

(51) 

<Pa(°a) 

(52) 
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To obtain a finite dissipated energy, it is necessary that 

dp 2 


(53) 


We study, however, in chapter III, correlation laws which do not 
verify this last inequality as is also the case for the law given by the 
equation (23) • We assume in these cases that the curve is approached in 
the region where the values of p are very small, that is, at the vertex 
of the correlation curve and that it must be rounded off there to make 
the first derivative zero. 


In order to determine directly a function suitable for representing 
the spectrum, one will pose, aside from ( 52 ), the conditions 


f <p A ( s ) ds ■ 1 
J 0 

(5M 

lim <?£(%) = § 

(55) 

ft— *0 


lim 9 a (Aa) = 0 

(56) 


ft— »«° 


The equations (12) to ( 15 ) which determine the relation between the 
spectrum and the correlation law will be written with the new notations 


OO 

<Pa(%) - § P cos(ft^s)R A (s)ds 


(57) 


p 00 

Ra(p) = / cos(sp)q> A (s)ds ( 58 ) 

J 0 

For a study of the longitudinal spectrum and correlation it suffices 
to replace p and A in the equations (48) to ( 58 ) by | and x; for 
study of spectrum and transverse correlation, these expressions are 
replaced by T] and y. 

In the study of the homogeneous and isotropic turbulence of a flow, 
one may apply the equation (39) which gives 
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The law of transverse correlation is given as a function of the law 
of longitudinal correlation by the relation corresponding to ( 36 ) which 
is written 

BjtCti) = ^(1 n) + \ (60) 

and to the equation ( 37 ) there corresponds 


H X (S) 


' til 


(s)ds 


(61) 


The equation (lt-2) gives 

_ d 2 g x (0) 

d * 2 


and 


1 1 ^5y(0) 

? " " 5 


( 62 ) 


When the turbulence is isotropic, one has = l/2ft x ; consequently. 


the equations (40) and (if-l) become 


"V( n y) ’ ^(“y) - n y 


^K) 

dfly 


(63) 


- 2 n x/ i/2 Q s *r i 48 


( 64 ) 


and represent the relation between the longitudinal spectrum and the 
transverse spectrum. 

The equation (^5) which gives the relation between the dispersions 
of the two spectra, will become 



x 


(65) 
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Thus, representing the spectra in dimensionless coordinates, one 
finds that in a flow of homogeneous and isotropic turbulence the dis- 
persion of the transverse spectrum equals half the dispersion of the 
longitudinal spectrum. Equation ( 65 ) represents the relation between 
the dispersions of the two spectra when the longitudinal spectrum is 
represented as a function of and the transverse spectrum as a func- 

tion of fiy. If one represents the two spectra as functions of the same 
variable, of & x , for instance, one will have the relation 



(65*) 


To the equation (43) there corresponds the very simple relation 


lx 



( 66 ) 


which shows that the dimension of the smallest eddies referred to the 
longitudinal- correlation length is equal to the inverse standard deviation 
of the longitudinal spectrum. 

If one studies the turbulence by following the particles in their 
motion, one determines the conditions which must be satisfied by the 
equations capable of representing the correlation law and the spectral 
function by replacing in the equations (48) to (56) p and A by 

and tL. The equations (21) and (22) of J. Kampe de Periet will be 
written 


^lKl) ' I f 0 

■a-M ■ So 


and will have the same form as the 


cos(a fcL s)RtL( s ) ds 

cos(sTL)q> tL (s)ds 
equations (57) and (58) • 


(67) 

( 68 ) 
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CHAPTER II 

METHODS USED FOR REPRESENTING THE EXPERIMENTAL CORRELATION 
CURVES OR SPECTRA BY ONE EQUATION 


7. Usefulness of the Representation of the Correlation 
Curve and of the Spectrum by One Equation 

Representation of the correlation curve by one simple equation may 
be very useful. Thus it is easy - when the law of longitudinal correla- 
tion between the components of simultaneous velocities is known - to 
calculate the law of transverse correlation for a flow of homogeneous 
and isotropic turbulence with the equation (60) . Inversely, if the law 
of transverse correlation is given, one finds , with equation (6l) , the 
longitudinal correlation. The dimension 7 \ (or l) of the smallest 
eddies responsible for the dissipation of energy by turbulent viscosity 
r»nn he determined by (62) for all functions which have a second deriva- 
tive at the origin. The greatest service the correlation equation 
renders is in permitt ing the calculation of the spectral function with 
the equation (57) without graphical integrations which are very time- 
consuming and have little accuracy. 

One may use the equations representing the correlation curve for 
study of the diffusion phenomena and for calculation of the measurement 
correction with a hot wire of nonnegligible length as well as for the 
correction of measurements with uncompensated hot wires. Moreover, it 
is desirable to be able to represent the turbulence spectra or the corre- 
lation curves by equations of the same general form the coefficients of 
which permit a comparison between different turbulent flows. 

If the spectrum is represented by an equation, the correlation law 
can be calculated with the equation (58). This correlation law may then 
be employed for other calculations. 


8 . Apparent Correlation Length 

When one has a certain number of experimental points at disposal 
and wants to determine the correlation law most convenient for repre- 
sentation of the test results, one will if possible begin with the cal- 
culation of the correlation length by planimetric measurement of the area 
bounded by the curve which best represents these points. When the corre- 
lation curve is such that negative correlation coefficients exist, that 
is to say, when it intersects the r-axis, or when this form is considered 
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possible, it will then he necessary to measure the coefficients R A (r) 
up to values of r sufficiently large to allow determination of 

with sufficient accuracy. It will frequently he viseful to employ, for 
the study of such correlation laws, what is called the "apparent" cor- 
relation length ... 


(ap) 

A 



R A (s)ds 


where Tq represents the smallest value of r for which the correlation 
curve intersects the axis of r (fig. 2). 


Using the dimensionless coefficients, one will have 

(ap) 

L A pPQ 

5i (s > d3 


(69) 


where pQ represents the smallest value of p for which R (p) becomes 
zero. Putting " A 


(ap) 


(ap) 


(ap) 


( a P) 


_ ( a p) _ r 


one will have 


,(ap) 


(ap) 


= a . n 


1 


(ap) 


X, 


(TO) 


For all correlation laws which do not give negative correlation 
coefficients, one will have 


(ap) 

L = L a and X. = 1 
A A A 


9 . Representation of the Experimental Results by One Equation 

(a) In the following chapters we give a large number of curves which 
represent correlation laws of different forms. For determining the most 
convenient one, one will begin by constructing the experimental curve 
representing g as a function of g, rj, or t, according to the case; 
we will draw it to the same scale as the curves given in the present report. 
Superimposing on these curves the experimental curve (drawn on tracing 
paper), one will see quickly which one is the equation which may best 
represent the experimental results. 
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After having chosen the form of the equation, one can employ t hr ee 
methods for calculating its coefficients. 

(a) By interpolation or by trial and error 

( 0 ) By taking as a basis the area bounded by the experimental curve, 
and the moments of different order of this area 

(7) When the correlation curve is represented by a polynomial, one 
can determine the coefficients by fixing beforehand the number of points 
of the experimental curve through -which one will have the calculated curve 
pass. One will then determine the coefficients by calc ula ting the coef- 
ficients of the polynomial. 

In order to apply the method* (b), one has to calc ulat e a certain 
number of moments of the area bounded by the correlation curve . Assume 


L (K) = f s K R(s)ds 

A J 0 


(71) 


to be the moment of the order k. One will have, in particular for the 
area bounded by the curve 


- (0) = / ?A (s)dS = 1 
A Jq A 


the moment of the first order will be 

( 1 ) 


L 1 = /'“ sR A (s)ds 

" A J 0 


and the moment of the second order 

( 2 ) 


"a = So s2 ^ s ^ ds 


After having selected the law which represents R A (p), one finds 

by integration the equations which give L^. On the other hand, one 

“A 

determines graphically or with a planimeter5 moments for the experimental 

^Certain planimeters (for instance, the apparatus of Koradi) measure in 
a single operation the area, the static moment, and the moment of inertia. 

Using such a device, one finds s imult aneously , and . 

"a "a "a 
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curve. Thus one obtains (K + 1) equations which permit calculation of 
the coefficients of the correlation law. Since the number of moments 
vised is never infinite and since the correlation curve resulting from 
this calculation can frequently be different from the experimental curve, 

K 

both giving the same values of L^, it is always necessary to verify the 
result . 

When the correlation curve intersects the axis and one does not have 
at disposal measuring results up to sufficiently large values of p, it 
may be useful to employ moments of the "apparent" surface of the correla- 
tion curve 



( 72 ) 


one will have, in particular 



(b) Measurement of the spectrum furnishes an experimental curve f (o>) 
as a function of on If the correlation length is known, the function cp(fl) 
is easily determined. In the opposite case one can calculate it if one 
assumes the form of the correlation law. For this purpose one traces f(oi) 
as a function of co using logarit hmi c coordinates, and one superimposes 
this curve on the curves corresponding to the selected correlation law, 
with consideration of the relation cuf(co) = £kp(&) . Thus, one finds the 
relation between .co and 0 and hence the correlation length, and one 
obtains at the same time the coefficients of the correlation law, and con- 
sequently also the equation of the spectral curve. 


The methods employed for representing the correlation curve by an 
equation may serve for finding directly the equation .which can represent 
the spectrum. In this case one will have for the moment of the order K 
the expression 


K 

F = 


s^CsJds 


and in particular 

F^ 0 ^ = f cp(s)ds = 1 



0 


( 73 ) 
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10. Experimental Results 

In table III we give a list of the experimental results on measure- 
ment of the correlation coefficients and of the turbulence spectra ■which 
will serve us in what follows for comparison with the correlation laws and 
with the spectral functions. Every experimental curve has a name (several 
letters followed by a number), and we shall indicate it in the future 
simply by this name. The largest number of experiments furnishes the 
longitudinal- correlation coefficients R x and transverse-correlation 
coefficients Ry. The curves KD.la, KD.lb give the correlation coef- 
ficients R^ and R^ between the components of the turbulent velocities 

perpendicular to the direction of the mean velocity. The experiment NPL.2 

v 

gives the correlation curve R taken from Lagrange's viewpoint, that is, 

■fc Ij 

considering the transverse turbulent velocities of the same particle at 
two different instants. This correlation curve has not been obtained by 
direct measurements but by performing a calculation starting from the 
turbulent -diffusion tests (ref. 6, p. 473). The experiment EGR.lc fur- 
nishes the correlation R m between the longitudinal turbulent velocities 
at a point which is displaced with the mean velocity of the flow. KD.lc 

V 

gives the correlation R between the transverse turbulent velocities at 

t 

a point fixed in space (with respect to the hydrodynamic center) . 

The spectra have been determined by making experiments at a point 
fixed in space, and concern the energy distribution of the longitudinal 
turbulence . Thus, one does not have at disposal either a longitudinal 
turbulence spectrum f x (o>) or a spectrum f^^cii) of J, Kampe" de Feriet. 

In table IV we give the values of the apparent correlation length 

and of the true correlation length L. The latter is given only 
for the experiments where the correlation curve intersects the abscissa 
axis. For the experiments NBS.2 to NBS.6, for which measurements have 
been made with positive and negative distances y, we calculated the 
correlation length by taking the average of the lengths measured for the 
two parts of the curve. In the same table we present also the moments 
of the first and of the second order of the "apparent" surface of the 
correlation curve. We did not calculate these moments for the experi- 
ments H.la, 2a, 3&> 4a for which their determination with sufficient 
accuracy would be difficult. 

„ r (o)i (ap) r (Di 

Figure 6 represents |_L V (J as a function of [L J for 

the correlation curves R x and Ry. On this figure only the experimen- 
tal points related to the wind-tunnel tests are given. It is found that 
these points can be represented rather satisfactorily by a straight line. 
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The points which correspond to the laws of the form B^(p) = exp(-|p|) 

and R^(p) = exp are placed rather close to that straight line. 

It is viseful to remember that in all tests made in the wind tunnels the 
turbulence is measured downstream of the grid. 


Assuming for the experiments for which one does not find a value of 

the apparent correlation length = L, one is able to give the ratio 

of the longitudinal correlation length and the transverse correlation 
length and to vferify the equation (39) • These ratios are given in table V. 


Iy 

For the experiments KD.la, lb we give the experimental value of which 

Xj 

X 

may be compared to the one given by equation (39 f )« 


In a flow of homogeneous and isotropic turbulence the values of these 
ratios must be equal to 2. When one studies the experimental results 
obtained downstream of a grid and demonstrates by the relation existing 
between longitudinal and transverse correlation that _ the turbulence Is 
homogeneous and isotropic, one has to understand that this property is 
approximate and concerns only the region in which the. tests have been 
made. One states, In fact, that the intensity of the longitudinal tur- 
bulence decreases when the distance from the grid increases. 
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CHAPTER III 

THE CORRELATION LAWS WHICH DERIVE PROM R A (p) = exp(- |p|) 


We study in this chapter the correlation laws of the form 


R a (p) - 


n=K 

A 0 + A n cos^cp) exp(- 
n=l 


c|p|) 


R A (p) = 


1 + 


A n . I n 

Z_ A n c |Pl 

n=l 


exp(-c|p|) 




(75) 


h a (p) 


n=K 


SI K exp(-Cn|Pl) 
n=l 


(76) 


dR A (°) 

For these laws, the condition — ^ — =0 is verified only in 

exceptional cases; the dissipation energy is therefore generally not 
finite. These laws are nevertheless of interest because they are rela- 
tively simple and can frequently represent satisfactorily the experimental 
correlation curves. Although it is not possible to apply them to the cal- 
culation of the value of l (or of A) , they can nevertheless serve for 
determining the spectrum of turbulence and for the application of Th. von 
Harman’ s law in a flow of homogeneous and isotropic turbulence. 


Since the aim of this study is to permit the selection of an uncom- 
plicated equation for representation of the experimental curve, we shall 
examine only the simplest examples of these laws. 


11. Law R A (p) = exp(-|p|) 

If one makes In the equations (7^) , (75) K = 0 and in equation ( 76 ) 
K = 1, one finds the simplest example of these laws 


R^p) * exp(-jpl) (77) 

The use of this equation for representing the transverse correlation 
curve has already been proposed by the National Bureau of Standards as we 
have recalled in the first chapter. H. L. Dryden also has utilized it. 
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by analogy, for representing the longitudinal correlation, and has deduced 
from it the spectral function (ref. 12 ). 

Making the Fourier transform according to (57) , one finds for the 
spectrum of the turbulence following a straight line in the direction A 
the expression , 

q&(OA) = | — (78) 

It -i . r, 


(a) Assume the longitudinal-correlation law corresponding to (77) 
be 

R x ( 6 ) = exp(-IU) _ (79) 

In a Flow of homogeneous and Isotropic turbulence, one may apply the 
equation ( 60 ) for determining the transverse-correlation curve which will 
be 


3yW = ‘ ^l T ll) ex I>(- | hi) 


and which one can represent as a function of applying ( 59 ) by 


RyU) = (i - £|M)"®(“IM) 


(80) 


( 80 ') 


The correlation curve R^.(t]) intersects the abscissa axis for the 

value t]q = 4 and becomes negative for the values larger than t), passing 

through the minimum: |Ry(T])l = -0.5 exp(-3) = -0.02^9 for tj = 6. 

J min 

Figure 7 represents the curves: R x (|), Ry(ri), and Ry(fe)* 


One can give the correlation law also with the apparent correlation 
length as a basis. One will Obtain X x = 1 and applying equation (69) 

one finds: X^ = 1 + exp (-2) = 1.1353* Taking (70) into account, one 

can trace By as a function of . The curve which represents the 

longitudinal spectrum corresponding to ( 78 ) 


cpx(ft) = | o (81) 

n 1 + <T 


is given in figure 8 . 



KACA TM 1436 


45 


(b) Assuming the transverse-correlation law 


By(Tj) = exp(-|T) |) ( 82 ) 

one finds, in a flow of homogeneous and isotropic turbulence, the 
longitudinal-correlation law 

§*(1) = ^[i - (1 + 2|!|) eX p(-2|i|)] (83) 

which can he expressed as a function of 7] by the relation 

5x(ri) = “ ( 1 + hl)exp(-hl)] (83-) 

n J 

The curves represented by the three last equations are given in fig- 
ure 9* To the .transverse-correlation law ( 82 ) there corresponds the 
transverse spectrum analogous to ( 78 ). In a flow of homogeneous and 
isotropic turbulence, one finds the longitudinal spectrum by applying 
equation (64) . One obtains 

= 1 1[| arc tan (!) - X ] w 

The two longitudinal spectra (8l) and (82) are compared in figure 8. 
Let us note that, by assuming the correlation law ( 82 ), one does not 
obtain the spectrum (8l) employed by H. L. Dryden (ref. 12), but the spec- 
trum represented by the equation (84) (-when the turbulence is homogeneous 
and isotropic). Besides, the two spectra do not differ greatly. 

(c) Comparing the transverse-correlation law (82) with the results 
of the experiment KBS.l, one finds (fig. 10) that the experimental points 
are relatively represented by the theoretical curve. A similar comparison 
made for the experiment NBS.7 gives a still better result (fig. 11). With 
the transverse-correlation law known, one can give the longitudinal cor- 
relation, for homogeneous and isotropic turbulence; it will here be deter- 
mined by the relation ( 83 ). The theoretical curve and the experimental 
points (KBS. 7a) have been drawn in the figure, for the longitudinal cor- 
relation R x (tj), and since the points are not placed on the curve, one 

can say imme diately that the turbulence of the flow in which the experi- 
ments KBS . 7 have been made is not homogeneous and isotropic. 
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12. Law 



n=K 

Aq + XZ An cos(mn cp) 
n=l 


exp(-c|p|) 


The equation (74) may he represented In a slightly different form 
by replacing 1 % which has an arbitrary value by integers n. One then 
obtains the law 


S A (p) = 



cos(ncp) 


) 

— 


exp(-c|p|) 


( 85 ) 


which comprises a Fourier series'. An application of this series, will, 
besides, be difficult because it will require, in general, employment 
of a rather large number of terms. We study here only the simplest laws 
•which have the form of the equation (74), and in the first place, the 
one which one obtains when Aq = 0 and K = 1. 

(a) Law I^(p) = exp(-c |p| )cos(mcp) .- To the correlation law of the 

form 


R A (p) = exp(-c|p|)cos(mcp) 


( 86 ) 


there corresponds the spectral function 


q&(r&) 

Integrating the 


.sT 1 + k . .1 

|c 2 + (me + fl) 2 c 2 + (me - ft) 2 J 
equation (86) from zero to infinity, one finds 


( 87 ) 


The spectral function is always positive. Since the equations (86) 
and ( 87 ) depend only on the absolute value of m, one can admit that m 
is always positive. - 

The correlation curves (fig. .12) have the form of a damped sine, 
the damping of which decreases with increasing m. When one has obtained 
an experimental correlation curve of this type, one might have diff iculties 
in measuring its correlation length with sufficient precision and then to 
represent it as a function of p, in order“to compare the experimental 
curve with the theoretical curves. It will be much easier to make the 
comparison on the basis of the apparent correlation length. The value Pq 
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for which the correlation curves intersect the p-axis is given by 

p n = — . Applying equation ( 69 ) , one finds the ratio of the apparent 
u 2mc 

correlation length and the true correlation length "which depends only on 
the value of the coefficient m 


X A = m exp 




and which is represented in figure 15 . Figure 13 gives the correlation 
curves as functions of p( a P). The abscissas of the intersection points 
of these curves with the axis of the p( a P) are represented as functions 

of m in figure 15 . Figure 16 gives the abscissa [p( a P)] R min of the 

minimum of the correlation coefficient as a function of m, likewise the 
value of this minimum 


k (ap) ] 


Emin 


X P Bmin 


= i — arc tani 
X me 1 


(-i + ”) 




m^ + 1 


• ^ arc tan 


(-i + ”)] 


The spectral cumres are represented in figure These curves pass 

through maxima, the values of which increase with m. Figure IT gives 


the value of 


to which corresponds the maximum of the spectral 


function and which one calculates by applying the relation 


- (m^ + l] 


c Pmax 


and also the value of 




(m 2 + l)' 


which one obtains by applying this 


relation to ( 87 ) . Since is always positive, the maximum of the 

spectral function corresponds to & A = 0, for all values of m smaller 
than \/l/3 . 

As has been said in the second chapter, the coefficients of the 
equation representing the correlation curve may be found by calc ulat ing 
the moments of different orders for the experimental curve and by comparing 
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them with the moments of the theoretical curve. In the actual case only 
one coefficient Is to be determined. Also, it is sufficient to know the 
moment of the first order. When one represents the correlation law as a 

function of p, one finds L ^ = 1 - m 2 , and by giving the curve as a 

function of p(ap) one obtains, applying equation (72) 



These two moments are given as functions of the coefficient m in 
figure 18 . 

(a) When one has, in a flow ^ of homogeneous and isotropic turbulence, 
a longitudinal-correlation law of the form ;7 

R^(g) = exp(-c|||) cos(mc|) (88) 


where 


one finds for the transverse correlation, applying the equation (60) : 

By(n) - [j ^ 1 " | c hl) cos (| met]) - £ mc|Ti|sin(! mc| rj (Jjexp | c|tj|)(89) 

One can express this correlation as a function of | , by applying 
to this last equation the relation (59) • 

( 0 ) When the transverse-correlation law in a flow of homogeneous and 
isotropic turbulence is .t 


Ry(T)) = exp(-c|T)l )cos(mcT]) 

one finds with equation (6l) for the longitudinal correlation 


(90) 




jm(||| + l)sin(2 mcl!l) - 

i(2||| + 1 - m 2 ) cos ( 2mc| )j exp ( -2c || | ) + i(l - m2) 


( 91 ) 
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( 7 ) Since the points of the experiment NBS .8 are placed along a curve 
the shape of ■which recalls the curves represented in figure 13 , we try to 
represent them by such a curve. For obtaining the coefficient m, one 

takes k (i >r p) in table IV and one finds in figure 18 m = 0.45« 
Transferring this value into equation (90), we find the correlation curve 

which is given as a function of p( a P) in figure 19 and which represents 
the experimental points quite satisfactorily. 

(b) law R A (p) = [Ap + A 1 cos(mcp)] exp(-c )p | ) .- When one makes in 
equation (74) K = 1, one obtains the correlation law of the form 


R A (p.) = [Aq + A-^ cos(mcpj] exp(-c|p|) 


(92) 


The spectral function which is calculated with the equation (57) 
will be 


<Pa(»a) 1 


2A 0 


Al 


c + ft c^ + (ft + mc)^ c^ + (ft - mc)^ 


( 93 ) 


For determining the coefficients of these equations, one must apply 
the conditions (48) and (55), and one obtains 


A 


1 



and 


A^m^ + 1 
m^ + 1 


which permits to give ( 92 ) and ( 93 ) as functions of only two coefficients. 
The coefficients Aq and m cannot be arbitrary, and their values are 
fixed by the inequation 


1 (5m.2 + 8) + \/l6m^' + 8Cta.R + 64 

- m 2 < A 0 ^ 9^2 


The ratio of the apparent correlation length and the true correla- 
tion length is calculated with the equation ( 69 ) * and one finds the 
relation 


X = 1 - 


Aq m 2 + 1 


mr~ 4- in 


Ap - 1 

Jfa - 


\jl - 2 Aq exp(-c|p 0 |) 


where 
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As a consequence , the correlation curve intersects the p-axis only 
when Aq < 0.5'. ... -- 

The correlation curve is monotonous in’ the interval 
(m 2 + l) - j/m 2 + 1 ^ A ^ (m 2 + 1 ) + ]/m 2 + 1 


and outside of this interval the abscissas of the ma.v-imn. and the rninlmn. 
are given by 

|PdR I = — 0X0 sin/ ' ■ — ■ . — 1 \ - arc sin ( — 1 — M 

I if =0 | mc |_ \Ao “ 1 Vm 2 + lj \jm 2 + iyj 

We give in figure 20 the value of Aq as a function of m verifying 
the conditions (-h-8) to (52). The admissible values of the coefficients Aq 
and m are divided into three types: 

dB A 

First, the coefficients for which — = ,< 0, that is. for which the 

dp v 7 3 

correlation curve does not present either minimum or maximum, with the 
curve descending continuously from 1 to 0 


Second, the coefficients for which 0 $ R. <: 1 but for which 

^ dp 

can become negative so that the correlation then presents positive maxima 
and minima 


Third, the coefficients Aq, m for which -1 ^ R^ ^ 1 which gives 
correlation coefficients presenting positive maxima and negative minima. 

Figure 21 represents the correlation laws as functions of p( a P) 
for a few values of Aq and of m, and figure 22 gives the corresponding 
spectral functions. One finds again in these figures _the curves already 
given previously since one obtains the law (77) when Aq = 1, and has the 

equation (86) for Aq = 0. 

These curves permit representing the experimental results by super- 
imposing the theoretical curves on points traced hs functions of p(ap). 
For trying to determine the coefficients Aq and m, using the surface 
moments of the curve, it would be necessary to give the moments of the 
first and of the second order. 

(a) When. the longitudinal-correlation law is of the form 


= 1^0 + A 1 cos(mc| i | )J exp(-c| | |) 
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one finds for the transverse correlation, in a flow of homogeneous and 
isotropic turbulence, the expression 

SyM - [a 0 (i - J cm) + A x (l - £ c|Ti|)cos(| men) - 

i A ± me |t] | sin^i mcl^ljjexp (- 1 c [ tj (95) 

which can be given as a function £ by 

By(!) = ^1 - c|li^R x (l) - \ A x mcl||exp(-c|i|)sin(mc|il) (95') 
(p) To the transverse-correlation law 


ly(tj) = (Aq + cos (merj )J exp ( -me | ^ | ) 

there corresponds the longitudinal correlation 


(96) 


R x (6) = 4 A l ffi2 + 4 Wmlil + l)sin(2mc || | ) - 
2s ^ 

i(2||| + 1 - m 2 )cos(2mc)||l - Aq^ + J| 


lexp(-2c|| | ) 


(97) 


which one can express as a function of tj by applying the relation tj = 2 | . 

( 7 ) Superimposing on the experimental points of the experiment NBS.l 
the curves of the figure 21 , one finds that one obtains a good representa- 
tion when the coefficient m = 0.5 and when Aq lies between 2 and. 4. 

We determined by trial and error that the best result corresponded to 
A 0 = 2.75* Figure 25 shows that the theoretical curve represents the 

experimental results for this test perfectly. 

Figure 24 gives the points of the experiment NBS.5 represented by 
the curve m = 0.5* Aq = 2. 

In figure 25 we represent the experiment H.4a by the longitudinal- 
correlation curve m 1.5# Aq = 0.6. This curve presents maxima and 
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minima for large values of | as could be predicted on the basis of fig- 
ure 20. Equation (95') permit's calculation of the corresponding transverse 
correlation curve -which represents very well the points of the experiment 
H.Ab for the values of | smaller than 1.2. One may therefore as stone 
that the flow for which the experiments H.4 have been made approximates 
rather closely a flow of homogeneous and isotropic turbulence . 


15. Law R^(p) = 


1 + 


n=K 1 

21 An c n (lPl) n 

n=l 


exp(-c |p I ) . 


(a) Law.- Rft(p) = [l + A 1 c Jp |] exp(-c |p| ) 

When one makes in equation (75) K ■ 1, one finds the correlation 
law of the form . ■. =_ 

+ A 1 c|pfjexp(-c|p|) (98) 

The corresponding spectral function is 






c 2 - 

—s + A 1 

Q (c 2 


- ft 2 " 

+ &f_ 


(99) 


By integrating the equation (98) from zero to infinity, one finds 


c = 1 + 

and prescribing the conditions (48) to (52), one calculates that the coef- 
ficient can lie between the limits 


■1 < ^ 1 


The abscissa of the intersection point of the correlation curve and 
the p-axis is given by 

,P01 “"Aifl + Ai) 

This shows that the curve intersects the axis only for < 0. One 

finds finally for the ratio of the apparent correlation length and of the 
true correlation length the expression 


1 


1 + Ai 
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The figures 26 and 27 represent the correlation law and the spectral 
function f.or several values of the coefficient A^_. When one uses coef- 
ficients of a value very close to - 1 ^ one must write them with a suffi- 
cient number of decimals because then the value of Aj_ has a very great 
effect on the shape of the curves. 

The minimum correlation coefficient is given by the relation 


L^ (p) ] inin * A i ^(tt) 

and corresponds to the abscissa 


P?min 


a i - 1 

Ai(l + A X ) 


These equations are valid only for Aj_ < 0. The value KL 

determines the. position of the maximum of the spectral function 
calculated with the use of the equation 


which 
may "be 


,c Rmax ^ 
which permits determination of 


1 + 3Ai , n 2 

ryl 1 + A l) 


max 


hy application of equation ( 99 )* 


This maximum lies on the cp-axis for all values of At > - 1. 

- 1 - 3 


In order to calculate the coefficient Ai by the method of moments, 
it suffices to know the moment of the first order which one finds by means 
of the equation ( 71 ) 


L (l) = 1 + 2A X 

" A (1 + Ai)* 

This method will be applicable only when < 1. The moment of first 

A 

order for the positive surface is obtained with equation ( 72 ) and is equal 
to 




(1 + A X )‘ 


t 


1 + 2A-i 


+ 
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The figures 28 to 31 represent the different characteristics of the 
correlation curves (98) and of the spectra corresponding to them. 

Let the longitudinal-correlation law be of the form 

5 X U) = (l + Ai cf|| )exp(-c||| ) (100) 

to which there corresponds in a flow of homogeneous and isotropic turbu- 
lence the transverse-correlation law 

Ry(ri) - |l + J(3A X - 1) c f T] | - ^ A x c^T|^j exp (- J c|tj|^ ( 101) 

which has the form of the laws given by equation (75). 

To the transverse-correlation law 


SyU) = (1 + Ai c|T)|)exp(-c|T}|) (102) 

there corresponds in a flow of homogeneous and isotropic turbulence the 
longitudinal correlation 

" x(l) = 2^2 j 1 + ^1 " [( X + ^ l)^ 1 + £c ^l > + ^ A 1 c 2 4^]exp(-2c||| )| 

(103) 


In order to represent the points of the experiment NBS.9 (fig. 32) 
by a law of this form, one may first try to apply the method of moments. 

(ap) 

= 0.847, given in table IV, as a basis, 

A^ in figure 31* For obtaining the moment of 


Taking the value of 
one finds the coefficient 




first order having this value, one may take either Aj_ = 0.6 or Aj_ = -0.5 . 

Superimposing on the points the figure 26, one finds that for the two coef- 
ficients the curves represent the experimental points gather satisfactorily. 
One sees, however, that a better result is obtained with the coefficient 
Ai = -0.25 (fig. 32). Besides, it will happen quite frequently that one 

obtains better results by trial and error than by the method of moments. 


For the experiment H. 4b, Li' T =0.762, hence, figure 31 gives 
Aj_ a* 0.95* One "sees that the transverse correlation curve corresponding 

to this coefficient represents very well the experimental points (fig. 33) . 
The longitudinal- correlation curve given by the equation (103) likewise 
represents very well the results of the experiment H.4a. As a result, the 
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flow for which these experiments have been made may be considered as 
homogeneous and isotropic. Figure 33 may "be compared with figure 25 
which also ■ represents the experiments H. 1 )-. In the first figure we repre- 
sented the measurements of the transverse correlation by an equation and 
deduced from it the longitudinal correlation by assuming a flow of homog- 
eneous and Isotropic turbulence. Regarding the second, we selected first 
the equation which represents the longitudinal-correlation curve and per- 
formed then the calculation for the transverse correlation. 

(b) Law.- R^(p) = [l + A x c |p[ + A 2 c 2 p 2 ] exp(-c |p| ) . 

By mftk-1 ng in equation (75) K = 2, one obtains the correlation law 
of the form 


R A (p) = [l + A x c|p| + A 2 c 2 p^]exp(-cfp| ) 


(104) 


The spectral function which corresponds to this correlation law is 
written 


%( 8 a) - f -rha 


c + 


c 2 _ n2 c 2 (c 2 - 3fl‘ 

1 + A-l u +2A 2 


C + £lr 


£) 

( C 2 + fi 2) 2 


(105) 


Integrating equation (10^), one obtains 


c = 1 + + 2 A 2 

and applying the conditions (48) to ( 52 ) one finds the limits for the 
coefficients of these equations 

-(l + 2A 2 ) < A X ^ 1 -1 < A 2 <c — ~ ^ ^ 

The abscissa of the intersection point of the correlation curve with 
the p-axis will be 

ill -Ai + jjj| - 4 a 2 

' P 0 1 c 2A 2 


and the ratio of the apparent correlation length and the true correlation 
length Is equal to 


= 1 - jj. + (A-l + 2A 2 )p 0 + A 2 cp^Jexp(-lc|p Q ) 
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Making the derivative of (10*0 equal to zero, one finds the abscissa 
of the point for -which the correlation coefficient is minimum or mavirmim 

1 2A 2 - A x t |/f% ■ 2A 2 ) 2 - 4A 2 (l - A X ) 

C OA. 


and which gives with the equation (104) the value of this coefficient. 

The above relationships permit tracing of figure in which we have 
indicated the limits within which the coefficients A^, A 2 may vary with- 
out causing the equation (104) to cease representing a correlation law. 

We have, moreover, defined the domain in which one finds monotonous cor- 
relation curves, furthermore the domain where the correlation curve pre- 
sents a negative minimum and fiAally the values of the coefficients for 
which the curve presents a positive maximum (without counting R^( 0 ) = 1 

which is not a maximum) . 

The figures 35 and 3 6 give the correlation curves and the spectral 
functions for -several values of the coefficients Aj_ and A 2 . One finds 

in these figures the curves presented already previously for the equa- 
tion ( 98 ) which is a particular case of the correlation law now being 
studied. 

When one knows the value of the moments of the first and second 
order, one can calcul a te the coefficients by solving the equations 

l/ - 6 L^ c 2 + 6 c - 2 = 0 ^ 


A 1 = -[^ (l)c2 “ 3c + 2 J A 2 '= |(c - A 1 - l) 

We have stated that the method of trial. and error gives very good 
results and is much faster than the method of moments. 

When in a flow of homogeneous and isotropic turbulence the longitudinal- 
correlation law is 

5 X (|) = JjL + A 1 c|i| + A 2 c 2 ! 2 ] exp(-c ||| ) (106) 

one finds, with application of von Harman 1 s law, for the transverse cor- 
relation the expression 

+ J(3 A 1 “ l)chl + |(^A 2 - A^c 2 !] 2 - ^ A 2 c^|ri| 5 Jexp(-| c |tj |) 

( 107 ) 
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This equation likewise lias the form of the law (75) , "but with K = 3« 
When, inversely, the transverse-correlation law is of the form 

Ry(Tl) = [l + A x c|ti| + A 2 c 2 |ii 2 |jexp(-clTi|) (108) 

one obtains in a flow of homogeneous and isotropic turbulence the longi- 
tudinal correlation 

R x (s) - 2 2^2 it 1 + 2*1 + Q- * U + 2c) III] exp(-2o||| ) - 

t( A l + 3A2)c 2 i 2 + 8A 2 C^ |l^ ljexp(-2c| II )^> (109) 

The figures 37 to b-2 show the results of several experiments and the 
curves which represent them. These curves have been obtained by trial and 
error. The results of the experiment KBS. 2 are represented by the 
transverse correlation law (108) with coefficients Aj_ = 0 and A 2 = 0.1. 
As shown in figure the curve intersects the abscissa axis passing 
through a minimum. Consequently, we have drawn this curve as a function 

of tj( a P) (fig. 37). 

The experiment KD.lc (fig. 38) concerns the correlation between the 
transverse velocities at the fixed point, and the experimental points may 
be represented by a law analogous to (104) , of the form 


?£( T l) = [l + A x c| 
or 

T 1 1 + A 2 c ^ t ‘ 
^ 00 

-Jexp(-c It|) 

* 1 ^ 

11 

H 

and 

O 

II 

R V ( s ) ds 
P 

with the coefficients Aj_ = 

-0.5, 

A 2 = 0 . 3 • 

One sees in figure 3^ that 


the correlation curve does not intersect the abscissa axis and that it is 
monotonous . 

The results of the experiment NBS.l have already been represented 
very well by the law ( 96 ), in figure 23* One can represent them as well 
by the correlation law (108), as shown in figure 39* 



58 


MCA TM 1456 


The measurements of the transverse correlation EGR.lb are represented 
in figure 40, "by the curve of the equation (108) with the coefficients 

and Ag = 0.5. This curve intersects the abscissa axis, and one represents 

it as a f unction of In a flow of homogeneous and isotropic turbu- 

lence, there corresponds to the transverse-correlation law (108) the 
longitudinal-correlation law (109) which one has also represented as a 

function of taking into account the relation 


(ap) = 2 5 


X 


The experimental points of the experiment^ EGR. la are not very well 
represented by this last law although some points are placed, on the 
curve R x . The experiments EGR.l have therefore been made in a flow, 

the turbulence of which is not fully homogeneous and isotropic but 
reasonably close to that state. 


The experiments KD.la and KD.lb give the . longitudinal and transverse 
correlations between the simultaneous transverse velocities. In figure -41 

we represented the correlation curve R V (experiment KD.lb) as a function 

_ y 

of T) v = — . If the flow in which these experiments were made were 


9 


homogeneous and isotropic 3 


£ = Sy 

would correspond to the equations (5) which would permit application of 

■y 

the equation (107) for calculation of the correlation curve R x . Since 

the experimental points KD.la deviate from this curve, the flow is not 
homogeneous and isotropic. 

Representing the results of the experiment H.2a by the longitudinal- 
correlation law (106) (fig. 42), one finds the transverse correlation for 
the homogeneous and isotropic turbulence , and one sees that the points of 
the experiment H.2b,do not correspond to the curve_ Ry. 

n=K . 

14. Law R A (p) = '2Z. ^ exp (-c n |p|J 

n=l ' ' 

The number of arbitrary constants one will have for the laws of 
this form is 2(K - 1). When K = 1, one has again equation (77) for 
which there is no arbitrary constant. We study here only the correlation 
law with K = 2. 



NACA TM 1456 


59 


(a) Law R A (p) = A exp(-cjp| ) + B exp(-3c| p ( ) . - Assuming in equa- 
tion ( 76 ) K = 2, A]_ = A, Ag = B, = c, and c 2 = pc, one obtains 
the correlation lav of the form 

R a (p) = A exp(-clpl) + B exp(-pc|p|) (llO) 
to which corresponds the spectral function 


- f- 


c^ + Q, c 


+ BP 


p I 2 c 2 + n 2 . 


(Ill) 


Integrating equation (llO) and applying the condition (50), one finds 


c = A + f B = 1 - A 

P 

which permits expression (110) and (ill) as functions of two constants A 
and p. 

The coefficient p may be larger or smaller them unity, but one will 
obtain exactly the same results by replacing this coefficient by its 
inverse. Assuming that 0 < p ^ 1, one finds for A the limits: 

1-0 1-0 

The abscissa of the intersection point of the correlation curve with 
the p-axis will be 

I P ° I " c(l - p) 102 (a - l) 

and since it must he positive, the curve intersects the axis only for 
A > 1. The ratio of the apparent-correlation length to the true- 
correlation length is 



The abscissa of the point for which the correlation coefficient is 
minimum, is equal to 


Bmin 


1 

c(l - P) 


log 




60 


NA.CA TM 1436 


and the value of this minimum is of the order of 

KL» - A 

The value of for ■which the spectrum presents a ma ximum may he 

calculated with the equation 

2 _ 0 2 (l - P 2 ) V-ABP - (a? 2 ♦ Bp) 

*»i A + SP 

and one obtains this maximum by substituting this value into equation (ill). 

Figure 43 gives the limits within which A may vary as functions of 
0. For the present law there are no correlation curves which present a 
maxi mu m as was the case for (104) . There are only two types of curves: 
those which intersect the axis and pass through a min imum value, 
those which do not intersect the axis and are monotonous. The correlation 
curves and the- spectral curves are traced in figures 44 and 45 for several 
values of 0 and of A. The different characteristics of these curves 
are given in figures 46 to 51. 




The moments of the two first orders are given by the relations 





When one knows their values for an experimental correlation curve, 
one can then determine the coefficients of equation (110) by solving the 
equations 


1 l ( 2 ) 

2 i 


i l (2) 
2 “ - 



0 = 1 - - c - , . • ■ A - 

c - C 2 l(D .... P - 1 - 

The moments of the "apparent" surface of the correlation curve are 
given by the equations 
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6l 


M 


(ap) 


= — Ja[ 2 - (a 2 + 2a + 2)b] + -^[2 - (a 2 p 2 + 2ap + 2)b P ' 

A 5 1 L - 1 p 5L 


5) 


or 


a = 


1-0 


log 


fcrh) 


and 


-H 1 ) 1 


i_ 

-P 


Figure 52 gives the moment of the first order 


m 


(ap) 


as a 


r rp^Cap) 

function of A for several values of p. Figure 53 represents |_L'- ; J 

as a function of [lU 3 (sp) . The curves corresponding to the different 
values of p coincide in this figure when A > 1. Thus, when A > 1 
and when the correlation curves intersect the abscissa axis, one can 
determine several correlation laws which give the same moments of the 
two first orders. 


Assume the longitudinal-correlation law to be of the form 


R x (g) = A exp(-c|g|) + B exp(-pc|g|) (112) 

In a flow of homogeneous and isotropic turbulence, there corresponds 
to this law the tr ans verse correlation given by the equation 

R^Xn) = A^l - i c|ii|^exp(-i c | T[| ^ + B^l - J pc| r\ |jexp Pc | n|^ 

( 113 ) 

When one has in a flow of homogeneous and isotropic turbulence the 
transverse -correlation law 


Ry(Ti) = A exp( -c |ti| ) + B exp(-Pc|ri| ) 
the longitudinal correlation is equal to 

IxU) = - \ 2 |a[i - (2c 1 1 1 + l)exp(-2c||| )J + 


(114) 


-^[l - (2pc||| + l)exp(-2pc[|| )] 


( 115 ) 
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One finds in table II, for the experiment NBS . 3 , the moments 


[id)]' ‘ 8P) =0.894 and ' l "} 


1.561. To these coordinates there 


corresponds in figure 53 the coefficient p = 0.25 with A < 1. Tracing 
then in figure 52 a horizontal, the ordinate of which is 0.894, one finds 
at the point of its intersection with the curve p = 0.25 a coefficient A 
which is approx ima tely 0.25 . The correlation curve corresponding to these 
coefficients actually represents the experimental points very satisfactorily 

(fig. 54). 


In the same manner one obtains for the experiment NPL.3 fo r ■which 

£l(i)J = 0.752 and = O.976, the coefficients p = 0.8 

and A = -3. The correlation curve which corresponds to these coeffi- 
cients is traced in figure 55 • 7 


This method is not successful, for the experiments H.la and NPL.la 
for which we shall determine the coefficients by trial and error. For 
the first time, we represent the points by the law (112) with the coeffi- 
cients p = 0.25 and A = 0.50 (fig. 5 6). Drawing the curve which 
corresponds to equation (113 ) > one sees that the points of the experi- 
ment H.lb are not to be found on this curve j this shows that the flow in 
which these tests have been performed is not homogeneous and isotropic. 


The results of the experiment NPL.la are represented by the 
longitudinal-correlation law (112) with the coefficients p = 0.25 
an d A = 0.25 (fig. 57)* The points of the experiment NPL.lb are placed 
very exactly on the .curve which corresponds to the- law (113)* Hence 
the result that the flow in which these tests have been performed is 
homogeneous and isotropic. 
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CHAPTER IV 


THE CORRELATION LAWS WHICH ARE DERIVED FRCM R A (p) = exp (-^ p 2 ) 


In this chapter we shall study the correlation laws of the form 


R a (p) = 


n=K 

A o + 21 ^ cos K °P) 
n=l 


exp(-c 2 p 2 ) 


(116) 


n=K 

? A (P) =1EZ A n exp (-c n 2 P 2 ) 
n=l 


(117) 


Since for these laws 


^°> - 0, 
a P 


it will be possible to determine 


by calculation the value of l (or of ~K) • As equation (66) shows, the 
relative dimension of the smallest eddies l x is equal to the inverse 
of the standard deviation of the longitudinal spectrum. 


15. Law R A (p) =exp(-g-p 2 ) 

When one makes in equation (ll6) K = 0 and in equation (117) 

K = 1, one obtains Gauss* curve 

R A (p) = exp (- ^ p 2 ) (118) 

Applying equation (57 ) , one finds for the spectrum the relation 

^(“a) -§«»(-&£) (119) 

which represents a curve of the same form as the correlation curve. This 
curve is given in figure 6l (for m = 0) . 

(a) When the longitudinal-correlation law is represented by the 
equation 


R x (i) = exp jj- | 2 j 


( 120 ) 
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one has, in a flow of homogeneous and isotropic turbulence, for the 
transverse correlation the expression 


- i 1 - & ’i 2 ) ei ‘»(-s "*) 


( 121 ) 


The curves representing these two equations are represented in fig- 
ure 58 . The curve which corresponds to (121) intersects the axis at the 
4 

point t]q = — whence: 

VJt 

of the point where the transverse correlation is minimum is equal to: 


X = erf(l) + -p= exp(-l) = 1,2581. The abscissa 
y VJt 


T Hni = Tn = 5*1915 and the value of this minimum is: 

R m j_ n = -exp (-2) = -0.1353* Applying equation. ( 62 ), one finds: Z x = 

(b) In a flow of homogeneous and isotropic turbulence where the 
transverse- correlation law is written 


RyU) = exp(- g. q 2 ) 

one finds for the longitudinal correlation the expression 


( 122 ) 


Bx(S) = -Kb. - exp(-rt|2) (123) 

U J 

The curves which correspond to these two equations are given in figure 59- 
For the dimension of the smallest eddies, the value ly = corresponds 

to this law. 

(c) The experiment EGR.lb can be relatively well represented by the 
curve (12), as figure 60 shows. The points of the experiment EGR.la are 
placed rather close to the curve which corresponds to the equation ( 123 ) • 
As a result, the flow rather approximates a flow of homogeneous and 
isotropic turbulence. This has already been seen when the test points 
of the experiment EGR.lb were represented by the equation (108) (see 
fig. 40) ; whereas it - had not been possible to calculate, for that last 
law, the value of ly algebraically, one obtains now immediately 
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A n 003 ( 1 % cp)J exp(-c 2 p 2 ) 

In studying these laws with K = 0 and K = 1, we did not find an 
application for the experiments given in table I. We present, neverthe- 
less, a few curves for K = 0 as well as the equations for K = 1. 

(a) Law R^(p) = exp(-c 2 p 2 )cos(mcp) .- When one 33iakes K = 0 in the 
the equation (ll6), one obtains the correlation law 

R^(p) = exp(-c 2 p 2 )cos(mcp) (124-) 


16. Law Ba(p) = Aq + 


= A 


n=l 


to which corresponds- the spectral function 


^ = - th & °) 


(125) 


Integrating the equation (124-) from zero to infinity, one obtains 


c = f exp (- 

Applying the conditions (48) to (53) , one finds that the coeffi- 
cient m may be arbitrary. One assumes, consequently, that m is always 
positive. 


The abscissa of the intersection point of the correlation curve and 
the p-axis is given by | p | = , and one finds the ratio of the apparent 

correlation length and the true correlation length by integrating (124-) 
from zero to Pq 


X A = 1 + 




— ^ 

(ni)I [2(n - npjjiVa 2 / 


In order to calculate this coefficient, it is easier to draw the correla- 
tion curves as f unctions of p and to determine L^ a P) by planimetering 
the area bounded by the curve up to its first point of intersection with 
the abscissa axis. 


The correlation curves and the spectral curves are represented in 
figures 6l and 62. 
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(b) Law R A (p) = [Aq + Aj_ cos (mcp )J exp (-c^p^. - To the correlation law 


g A (p) = |a 0 + Aj_ cos(mcp)J exp(-c 2 p 2 ) ( 126 ) 

corresponds the spectral function 

’’aPa) - VJJ_ A o + A i e ^(-T) ch fe 4 e5 ®(- &>) 

Integrating the equation (126) one calculates the value of the coef- 
ficient c which is 


<= = £J 


*o[f- - e3 ®(- Tr)] + e *p(- t) 


Applying the conditions (48) to (53) one finds that the coefficient Aq 
may vary within the limits 


0 < A 0 + 1 

m 


!The abscissa of the intersection point of the correlation curve and 
the p-axis is given by the relation 

Pol -E" 6 oos ( a^Tt) 

which shows that this curve intersects the axis when Aq < 

n=K - 

17. Law R^(p) = XI A n ex P(" c n a2 ) 


We shall study for these laws the particular case where K = 2 with 
the assumptions Aj_ = A, A<^ - B, c-^ = c, Cq = Pc. 


law 


(a) Law R A (p) = A exp(-c 2 p^) + B exp (~p 2 c^p^) To the correlation 
R A (p) = A exp^-c 2 p 2 ) + B exp(-p 2 c 2 p 2 ^ (128) 
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corresponds the spectrum 

*a(«a) - 

Integrating the equation (128 ) , one finds 


A exp (-&) + 1 


2 „ 2 , 


(129) 


c " “ 1 ) +1 ] B = 1 - A 

and applying the conditions (48) to (53 ) > one obtains the limits within 
which the coefficient A may vary 


0 < A ^ 


1 


1-0 


The different characteristics of the correlation curves are given 
by the relationships 



X A “ S A erf (epo) + erf (pcpo)J 



The figure 63 shows the values of the coefficients A and p which one 
can assume for the correlation law (128) • The curve intersects the axis 
when A > 1 while passing through a minimum value; it is positive and 
monotonous for 0 < A < 1* 

Figures 6k- and 65 give the correlation curves and the spectral 
curves for some values of A and of p. Figures 66 to 69 give the char- 
acteristics of the correlation curves, referred to the apparent correla- 
tion length. 
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The moments of the true area bounded by the correlation curve are 



and for the "apparent" surface one has 



The moment J_L V 'J is given as a function of A and of p in 

figure 70, |>) is given as a function of jj,(l)J i n fig_ 

ure 71 • For A > 1, one has one and the same curve for all values of 3. 

When the correlation curve does not cut the abscissa axis, one can 
also determine the coefficients by calculation, by solving the equations 



Let the longitudinal-correlation law be of the form 


= A exp(-c 2 £ 2 ) + B exp(-p 2 c 2 ! 2 ) 


( 130 ) 
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In a flow of homogeneous and Isotropic turbulence the transverse 
correlation which corresponds to this law will be 

Ry(n) = a(i - i c 2 Ti 2 )exp (-± c 2 Ti 2 j + b(i - ^ P 2 c 2 T) 2 )exp(- J P 2 c 2 ti 2 ) (131) 

and the relative length of the smallest eddies "which is calculated by 
application of the equation ( 62 ) is given by the relation 


*x*- 


2c 


: (A + Bp2) 


This length is represented as a function of A and of j3 in 
figure 72 . 

When the transverse-correlation law is 

Ry(T}) = A exp (-c 2 tj 2 ) + B exp(-P 2 c 2 Tj 2 ) 
one then finds in a flow of homogeneous and isotropic turbulence 

r 

H x (5) = A [ X - exp(.4c 2 £ 2 )]-^|l - exp(-4p 2 c 2 | 2 ) 

V- 

and equation ( 62 ) gives 


(132) 


(133) 


(13*0 


v 2 = 


c 2 (a + BP 2 ) 


(135) 


The relative length ly is given as a function of A and of p in 
figure 73* 

I" ,, H (ap) 

To the experiment H.2b correspond the moments I IA 'J 


= 0.893 




and = 1.410. The point drawn with its coordinates in fig- 

ure 71 is placed rather close to the curve p = 0.25. Then one traces a 
horizontal line with the ordinate 0.893 In figure 70 and finds that it 
intersects the curve p = 0.25 at two abscissa points A = 0.55 
and A = 0.88. In figure 71* the experimental point lies near to the part 

of the curve p = 0.25 closest to the axis Therefore one 

must in figure 70 take account of the intersection point with the branch 
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of the curve closest to the same axis. One thus obtains the two coeffi- 
cients p = 0 . 25 , A = 0 . 55 . The experiment is represented by a point 
which gives a negative correlation coefficient. In order to have a law 
which represents" a correlation curve which passes through negative 
values, it would be necessary that A > 1. So one sees that the point 
which corresponds to the moments of- the experimental curve is very far 
distant from the curve A > 1 (solidly drawn in fig. 71) • Thus it will 
not be possible to find coefficients for the equation ( 113 ) in such a 
manner that the curve represents the points of the experiment H.2b and 
intersects the abscissa axis. 

The transverse-correlation curve calculated with the coeffi- 
cients A = 0.55, P = 0.25 and the experimental points H.2b are drawn 
in figure 74. 


Comparing the points of the experiment H.2a with the curve R x cal- 
culated with the equation (134), one sees that the turbulence is not 
homogeneous and isotropic. This has, besides, already been seen in fig- 
ure 42. The relative dimension ly can be determined by the equation (135) 

or by the figure 73 as soon as one assumes the turbulence to be isotropic. 
In this case one would obtain ly = O. 632 . 


The moments corresponding to the experiment H.4b are [i/ 1 )] = 0.762 


and J^l(2)J _ 1 . 000 . One finds the coefficients p = 0.25 and A = 0.3 

in the same manner as for the preceding experiment. The transverse- 
correlation curve corresponding to these coefficients represents in effect 
satisfactorily the experiment - H. 4b, and the longitudinal-correlation curve 
calculated by means of the equation (134) overlaps with the points of the 
experiment H . 4 b; this shows that the turbulence of the flow in which these 
tests have been made may be considered homogeneous and isotropic. One 
finds as the ratio of the dimension of the smallest eddies and of the 
longitudinal-correlation length ly = 0 . 621 . 


Superimposing the points of the experiment KBS. 9 on the curves of 
the figure 64, one finds very quickly that they can be represented by the 
equation (133) with the coefficients p = 0.25, A = 0.5 (fig. 76 ) • One 


obtains easily, by calculation or on the basis of figure 73 , the value 

ly = 0.619. 


Likewise, the results of the experiment KPL.3 are represented by the 
equation (133) with the coefficients p = 0.5 and A = 0.5 (fig* 77) • 
One finds, for these coefficients, ly = 0.952. 


r 
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CHAPTER V 

THE CORRELATION LAWS DEVELOPED IN A SERIES 
OF HERMITE POLYNOMIALS 


Tn reference 5 > J. Kamp e de Feriet represents the correlation curve 
"by an equation which involves an expansion in a series of Hermite 
polynomials .4 With the notations we are -using in the present report, the 
correlation law is written 


E&(P) 


n=k 

i+z: 

n=l 


B n H 2 n 




( 136 ) 


where 


H^Cs) 


(-l) n 2 nl_ 1 . 2 n s 2 + 2 2 - D B ^ - 

gn nl 21 4i 


. (_i) n 2 an 


2 n 


(2n)T 


and for the first polynomials one will have 


HgCs) = s 2 - 1 


Hg(s) = s 6 - 15s 4 + 45s 2 - 15 


Hj^(s) = s^ - 6 s 2 + 3 


Hq(s) = s 8 - 28s 8 + 2108^ - 420s 2 + 105 


To the correlation law ( 136 ) corresponds the spectral function 

2 ^ 

2 


■ to 

which involves an expansion in a power series. 


n=k f ci 

1 + 21 (-i)bJ-^ 

n=l \2c 


exp 


(-a 2 ) 


(X37) 


18. Law 


n=K 

1 + 

n=l 


A n c 2n p 2n 


exp 


(-= 2 p 2 ) 


It would be much more convenient to employ an equation for repre- 
sentation of the correlation curves in which the development In a series 
of Hermite polynomials Is replaced by a power series 


4 The study of the Hermite polynomials has been carried out by 
J. Kampe de Feriet (ref. 18). 
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5 a <p) - 


n=k 

21 

n=l 


1 + > c^p 211 


exp 


(-A 2 ) 


(138) 


The relationship between the coefficients of these two equations has 
been given by J. Kanrpe de Feriet and is written 


m - k ~n 

1 v - 2 1 (n + 2ni )l 

An “ et n |ro 2^r — srr B n + 2n i 


V 


n-, 


k-n 


(139) 


®n = 

“l* 0 




Applying equations (57) to (138) , one finds the spectral function 


(”A "a) - + 2 <-U“ £ 3 ) 


(14-0) 


which comprises a series of Hermite polynomials. 

Integrating the equation (138) from zero to infinity, one obtains 


c 




n=l 


• (2n - 1) 

A = 


(141) 


The inequality (53) leads to the general condition for all laws of 
the form (138) 

A x < 1 

In the equation (14-0) the coefficients, the indices of which are 
even numbers, are preceded by a negative sign. In order to confirm the 
condition (52) , It Is necessary that Ak > 0 when K is even and that 
Aj£ < 0 when K Is odd. In the study of equation (138) one sees that in 
this last case the correlation curve R^(p) will Intersect the p-axis. 

For determination of the other conditions which must be set up in order 
to make the equation (138) really a correlation law, one must first know 
the value of K and then apply the equation (52). 


f 
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The simplest law of the form given by ( 15 8 ) is obtained when 
K = 0. One then has Gauss' curve which has already been studied, equa- 
tion (118). 

(a) When in a flow^of homogeneous and isotropic turbulence the 
longitudinal-correlation law is of the form 


5x(S) - 


n=k 


2 n. 2 n 


1 + > A n c | 


exp 




(142) 


one finds, with application of equation ( 60 ), the transverse-correlation 
law 

R y (Ti) = . 1 + J(2A X - 1)c 2 t! 2 + $1 ^[(n+ 1)^ A n _J c 2 ^ 2 * 1 . exp(-£ c 2 n 2 ) 
^ ^ (145) 

which is a law of a form similar to (l^ 2 ) but with other coefficients. 

The equation ( 62 ) gives for the dimension of the smallest eddies 
responsible for the energy dissipation by turbulent viscosity the 
relationship 


*x = 


J /2(l - A x ) 

(b) To the transverse-correlation law of the form 

n=k 


(144) 


RyU) “ 


1 + > 

n=l ; 


exp 


(.cV) 


( 1 ^ 5 ) 


corresponds in a flow of homogeneous and isotropic turbulence the longi- 
tudinal correlation 


- xit) " 4c 2 | 2 


and the equation ( 62 ) gives 


1 / o 2 \ ^ 

= — ±— 4l - exp (-4c | ) + > nIAn 


n=l 


1 - 


(4c 2 ! 2 ) 


n l 


t — exp 


(-4c 2 | 2 ) 
(146) 


ly = 7 1 ■■ (1^7) 

op. - A x 
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19. Law R a (p) = j^. + A 1 c 2 p^jexp(-c 2 p 2 ) 

When in equation ( 138 ) K = 1, one obtains the correlation law 

R^(p) = [l + A-l c 2 p 2 ]exp(-c 2 p 2 ) 
to which corresponds the spectral function 

*A) ■ ifcf + 1 - f A i S) 


and where 


-#[ 1 + l A iJ 


(ikQ) 


(149) 


We have already shown that the condition A^ < 1 must be verified 

whatever the value of K may be. On the other hand, in the actual case, 
with K being odd, it is necessary that A K - A i < 0. Finally, the 

inequality ( 52 ) , applied to the spectral function (l49) , gives A^ > -2 
whence there results the condition 


-2< A 1 ^: 0 

The correlation curve always intersects the p-axis. The various 
characteristics of this curve as well as those of the spectral curve are 
given by the relationships 




X A = erf ( cp o) “ 2 A 1 Po ex p(- c2 P 2 ) 
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The figures 'JQ and 79 represent the correlation curves and the 
spectral curves for some values of A, and the figures 80 to 83 give the 
values of the expressions calculated above. 

(a) To the longitudinal-correlation law 

Sx-Cl) = [l + A-l c 2 ! 2 ] exp(-c 2 | 2 ) (150) 

there corresponds in a flow of homogeneous and isotropic turbulence the 
transverse-correlation lav 


5y(tl) = [l + J(2A X - l)c 2 Ti^jejq? c 2 i) 2 ) (151) 

(b) When in a flow of homogeneous and isotropic turbulence, the 
transverse -correlation law is of the form 

Sy(n) = [l + Ai c 2 t]^J exp (-c 2 tj 2 ) (152) 

the longitudinal-correlation law is 


kc 


f 

J— 

• -t 

,1 + A X - [ 

1 + + ^A-l c 2 |^J exp (-c 2 | 2 ) . 


-/ 


(155) 


20. Law Ba(p) = j^- + A^_ c 2 p 2 + A 2 c^p^J exp ^-c 2 p 2 j 

When in the equation ( 138 ) K = 2, one obtains the correlation lav 
5 a (p) - [l + A 1 c 2 p 2 + A 2 c^p^J exp (-c 2 p 2 ) (154) 


to which corresponds the spectral function 
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. <i55) 

and where 

c = -f£ + i A i + i A 2] 

The different characteristics of these curves are given by the relationships 



a %^ - |( A 1 + 5A 2 ) * /( A l + 5A 2 ) 2 - A 2 (* + 6A 1 + ^ 2 )] 

Applying the conditions ( 26 ) to (34), one finds for the coeffi- 
cients A-^, A 2 the limits indicated in figure 84. .The figures 85 and 86 

give the correlation curves and the spectral curves, and the figures 87 
to 90 represent the values of some characteristics of these curves. 



NACA TM 1^36 


77 


CHAPTER VI 
TURBULENCE SPECTRA 


21. Equations Suitable for Representing a Turbulence Spectrum 

We have shown in the first two chapters what the conditions are which 
must be imposed on an equation so that it can represent a turbulence 
spectrum. In the following three chapters we have given the spectral func- 
tions which correspond to a certain number of correlation laws. One can 
evidently perform the opposite operation by first assuming the equation of 
the spectrum and then calculating the correlation law. The simplest equa- 
tions one can propose for representation of the spectra are 


cp = A expC-BG^) cp = — — - cp = A exp ( -B 1 2 1 ) 

1 + Bfl 2 

See equations ( 78 ) and. (118). Writing them for the G. I. Taylor 
spectra, one will. have 

<P 0 (flt) = | exp ^ Gt 2 ) (156) 


and 


‘ftW'l— Va < 157) 

1 + Sit 

Imposing the conditions ( 52 ) and (5^) to ( 56 ), one finds for the 

third 

<Pfc(°t) = ! exp(-f K|) (158) 

and the corresponding correlation lav is written 


5t( T ) = 




( 159 ) 


In figure 91 we compare these three spectral curves ^ and in fig- 
ure 92 the correlation curves which correspond to them. 
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22 . Experimental Spectra 

For representing an experimental . spectrum by an equation, one may 
employ the same methods as for the correlation laws. Notably, one may 
search for this equation by trial and error after having superimposed 
the experimental points on the spectral curves drawn in the n ume rous 
figures given in this report. 

When the correlation length is not known and the experiment, con- 
sequently, does not give the curve but only f t (o>), one plots, 

in this case, on logarithmic coordinates Uf (cn) as a function of 

and then superimposes it on the curves noting that 




Figure 93 represents the results of the experiment NPL.^a as well 
as the spectral curve corresponding to the equation (129) where A = 0.6 
and p = 0.15. It should be noted that, in representing the correlation 
curve ( 128 ) corresponding to this spectrum, one will obtain a curve which 
does not represent the points of the experiment NPL.lb*. Thus, it is pos- 
sible to show that the verification of the G. I. Taylor equations of 
which we have spoken in chapter I, section 3 is not exact. This indicates 
how difficult it Is to make verifications of this type. 
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CHAPTER VII 
CONCLUSIONS 


1. There is much confusion regarding the numerous correlation coef- 
ficients and spectra studied by the various authors. We hope to have 
shown the difference between these diverse factors by the use of appro- 
priate notations. These notations may seem relatively complicated but 

we are of the opinion that one must not shy away from a complicated nota- 
tion which may make the ideas clearer. 

2. By studying homogeneous turbulence for which, moreover, the time 
averages of the turbulent velocity are equal to the space averages, we 
show that the correlation coefficients between the simultaneous turbulent 
velocities at two points placed along a straight line of given direction 
and a given distance apart are the same whether the study is carried out 

in the Euler or in the Lagrange system. 

3. We introduce in this report what we call the longitudinal spectrum 
of turbulence which is obtained by harmonic analysis of the curve repre- 
senting the components of the simultaneous turbulent velocities u* along 
a straight line parallel to the direction of the mean velocity. The trans- 
verse spectrum is obtained by study of the simultaneous u* along a 
straight line perpendicular to the direction of the main velocity. 

4. In the study of the equations of G. I. Taylor which establish the 
correlation between the simultaneous u* and the spectrum measured at 
the fixed point (ref. 7 ), we show that there exists a coefficient the 
value of which may serve as a criterion for the legitimacy of employing 
these equations. This coefficient 

must be very small; only then the G. I. Taylor equations can be considered 
exact. 

5* We give several important equations derived from Karman's law for 
a flow of homogeneous and isotropic turbulence • 

6. The representation of the correlation laws and the spectra by 
empirical equations may be of very great service. We made a very large 
number of applications to the experimental results, especially regarding 
the correlation curve. When the equation which represents the longitudinal- 
correlation curve is given, one may easily obtain the transverse-correlation 
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curve "by applying Karmen* s law, and if one has the experimental points 
for the two curves, one can verify whether the flow in which these tests 
have been made is homogeneous and isotropic. On the other hand, one can 
rapidly calculate the turbulence spectrum by making a Fourier transform. 
By the same method, on the basis of the transverse-correlation law or of 
the spectrum, one can study other characteristics of the turbulence. 

7. The functions of the form 

I^(r) = exp(-K|r|)cp(r) 

are generally more satisfactory for representation of the correlation 
curves than the functions 

I^(r) = exp (-K^r^) cp ( r ) 

One succeeds in a very good representation of several results, employing 
only two arbitrary coefficients in these functions. 

8. The application of Hernite polynomials (ref. 5) to the representa- 
tion of experimental correlation curves seems to be rather difficult since 
it is necessary for the representation of an experimental curve to choose 
more than two arbitrary factors. 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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TABLE H 

COMPARATIVE TABLE FOR SYMBOIS RELATINO TO THE CORRELATION AND TO 
THE SPECTRA USED BY DIFFERENT AUTHORS 

































































TABLE m 





Distance from 




Homs 

Measuring 

result! 

Grid mesh 
In cm 

the mo souring 
point to the 
grid In mesh 

Mean velocity, 
cm/s 

Remarks 

References 




lengths 





Correlation curves j 

HBS.l 

Ry(y) 

12.70 

4o 




Ref. 11, fig. 21 

BBS. 2 

•V<y) 

.63 

4o 

1,220 



Ref. 11, fig. 5 

MBS, 3 

%(y) 

1.27 

40 

1,220 


One will use the 

Ref. 11, fig. 5 

nm.br 

%(y) 

2.54 

40 

1,220 


experimental 

Ref. 11, fig. 5 

UBS. 5 
BBS .6 

BjrW 

By(y) 

8.25 

12.70 

40 

40 

1,220 

1,220 


points without 
applying the 
correction taking 

Ref. 11, fig. 5 
Ref. 11, fig. 5 

HBS.Tft 

HBS.Tb 

%W 

By(y) 


2.5* 

40 

1,220 


into account the 
length of the 
hot wire 

Ref. 11, fig. 28 

BBS. 8 

%(y) 

.63 





Ref. 11, fig. 23 

BBS. 9 

Bjr(y) 

2.54 





Ref. 11, fig. 23 

BBS. 10 

Bjrto 

12.70 






H.la 

H.lb 

B.W 

.SjrW 

j 

1.27 

32 

610 


Ref. 9, fig. 17 

H.2a 

H.2b 

*iw 

Ry(y) 

1 

j 

1.27 

?r. 3 

610 


Ref. 9, «g. 17 

H.Ja 

M*> 

1 

.6} 

28 

610 


Ref. 9, fig. 12 

B.Jb 

iV(y) 





H.4a 

H.4b 

Ry(y) 

1 

j 

.65 

83 

610 


Ref. 9, fig* 12 

BFL.1& 

HFL.lb 

BxW 

Ry(y) 

l 

7.62 

27.3 



Curve calculabed by 

Ref . 7, fig. 4 

HPL.2 

RtL(b) 

2.29 

25.3 

610 


means of diffusion 
measurements 

Ref. 6, p. 474 

KPL. 3 

Sy(y) 

2.29 


763 


Ref. 6, fig. 1 

EDR.la 

H*W 

] 





Tests made in a 

Ref. 16, fig. 3 

Hffi.lb 

Ry(y) 

Br(h) 


.93 

24 



hydrodynamic 

EQR.lc 

J 




channel 


ED. la 

£(x) 






Tests mde in a 


ED. lb 

«J(y) 




19.8 


hydrodynamic 

channel 

Ref. 15, fig. 3, 5 

KD.lc 

Bj(h) 









Turbulence spectra 

BBS. 11 

f t (®) 

2.54 

40 



Ref. 12, fig. 1 

BBS. 12 

ft(®) 

2.34 

160 



Ref. 12, fig. 2 

BBS.4e 

f t (ro) 



*57 



BPL.4b 

ft(®) 



610 



NPL.4c 

ft(®) 

7.62 

27.5 

762 


Ref. 17, table II 

hflM 

ft<®) 



915 



WL.ko 

*t(“) 



1,067 



M.l 

f t(®) 



100 

u 

Experiments near 
the wall 

Ref, 13, fig. 2, 3 
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TABLE IV 





L (ai>) 

L 

[L(D] (aP) 

[l(2)] (aP) 

H. 2 a 

H. 2 b 

1.800 cm 

0.783 cm 

0.770 cm 

0.893 

1.410 

H. 3 a 

H.3b 

0.635 cm 
0 . 3^0 cm 

0.334 cm 

0.911 

1.527 

H.4a 

H.4b 

0.980 cm 
O .558 cm 


0.762 

1.000 

HFL.la 

2.74 cm 


1.125 

2.350 

HPL.lb 

1.55 cm 

1.43 cm 

0.967 

1.718 

HFL .2 

0.0044 a 


1.157 

2.687 

NPL-3 

0.443 cm 


0.752 

0.976 

KD.la 

2.34 cm 


0.986 

0.905 

KD.lb 

2.90 cm 


1.048 

1.068 

KD.lc 

0.208 s 


I .250 

2.895 

EGR.lb 

; 

4* 


O .678 

1.578 



L (ap) 

L 

[l«] M 

[L( 2 )] (al) 

KBS .1 

2.72 cm 


1.039 

1.907 

KBS. 2 

0.251 cm 

0.245 cm 

0.909 

1.242 

KBS .3 

0.520 cm 


0.894 

1.561 

KBS. 4 

0.795 cm 


0.794 

1.194 

KBS. 5 

1.995 cm 

1.945 cm 

0.952 

1-563 

NBS.6 

2.84-5 cm 


1.032 

1.889 

KBS. 7a 

0.945 cm 


0.960 

1.774 

KBS. 7b 

0.682 cm 


0.957 

1.598 

KBS. 8 

0.240 cm 


0.863 

1.634 

KBS. 9 

O .768 cm 


0.847 

! 1-514 

KBS. 10 

0.82 cm 


1.008 

1.832 

E.la 

1.358 cm 




E.lb 

O .634 cm 

O .634 cm 

0.928 

' 1.587 


y 

In reference 16 the author has not given the measuring units . 
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HkSLS V 


Exp. 

MBS . 7 

H.l 

E.2 

H .3 

H.li- 

NEL.l 

Wiy 

1-39 

2.11 

2 . 314 . 

1.90 

1.76 

1.92 


Exp. 

KB.1 

Jj/S 

1.2k 


A 
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Figure 5. 
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Figure 7. 
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Figure 13, 
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Figure 21. 
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Figure 28. 



Figure 30. 





Figure 31 
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Figure 35. 



Figure 36, 
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Figure 39. 





Figure 41, 


Figure 42, 
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Figure 44. 



Figure 45 
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Figure 59 
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Figure 64. 



Figure 65, 









Figure 66. 



Figure 68. 



Figure 70 
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Figure 80. 


Figure 81. 
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Figure 83 
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Figure 92. 
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Figure 93. 


NACA - Langley Field, V*. 



